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ABSTRACT 

A  theoretical  analysis  is  developed  treating  the  "Pump-Jet"  pro¬ 
pulsive  unit  comprised  of  stator,  rotor  and  enshrouding  nozzle  by  using 
the  unsteady  lifting  surface  theory.  The  analysis  takes  into  account 
accurate  geometry,  realistic  flow  conditions  and  hydrodynamic  interactions 
between  all  lifting  surfaces  of  finite  thickness  distribution.  The  system 
is  assumed  to  be  immersed  in  a  non-uniform  flow  of  an  incompressible  fluid. 

Expressions  have  been  developed  for  loadings  on  all  interacting 
surfaces  and  corresponding  resulting  forces  evaluated  at  proper  frequencies 
dictated  mainly  by  those  of  the  rotor. 
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INTRODUCTION 


Previous  investigations  at  Davidson  Laboratory  have  been  concerned 

with  the  adaptation  of  linearized  unsteady  lifting-surface  theory  to  the 

1  2* 

cases  of  a  marine  propeller  operating  in  a  nonuniform  inflow  field,  ’ 
of  counterrotating  propeller  systems, and  of  ducted  propel lers,'*’^ 
where  the  exact  geometry  of  the  systems,  the  realistic  inflow  conditions 
and  the  mutual  interaction  of  all  lifting  surfaces  are  taken  into  account. 

In  the  case  of  the  single  propeller  with  enshrouding  nozzle,  both 
accelerating  and  decelerating  ducts  were  discussed,  the  accelerating  (Kort) 
nozzle  offering  the  advantage  over  conventional  propellers  of  increasing 
the  flow  rate  through  the  propeller,  reducing  the  loading  and  thereby  in¬ 
creasing  the  efficiency,  and  the  decelerating  type  of  reducing  the  flow 
rate,  thus  delaying  cavitation  inception  and  lowering  noise  level. 

The  present  study  treats  the  pump-jet  configuration,  which  is  a 
type  comprised  of  stator,  rotor  and  enshrouding  nozzle.  The  stator  vanes, 
in  addition  to  their  structural  support  of  the  nozzle,  are  presumed  to 
homogenize  the  inflow  to  the  rotor  blades,  reducing  further  the  vibratory 
loading  and  resulting  forces  and  the  radiated  noise.  To  assess  the  advan¬ 
tages  or  disadvantages  of  the  system,  a  theoretical  analysis  and  correspond¬ 
ing  computer  program  are  developed  which  will  reveal  the  steady  state  and 
vibratory  characteristics  of  this  propulsive  device  as  a  function  of  various 
geometric  parameters  of  the  system. 

This  study  was  co-sponsored  by  the  Naval  Sea  Systems  Command  Explora¬ 
tory  Development  Program  and  General  Hydromechanics  Research  Program  under 
Contract  N00014-77-C-0298,  administered  by  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center. 


‘Superior  numbers  in  text  matter  refer  to  similarly  numbered  references 
listed  at  the  end  of  this  technical  report. 
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STATEMENT  OF  THE  INTERACTION  PROBLEM 


A  pump-jet  configuration  comprised  of  stator,  rotor  and  enshrouding 
nozzle  is  immersed  in  a  nonuniform  flow  of  an  ideal  incompressible  fluid. 
Figure  1  shows  the  relative  location  of  each  member  and  the  corresponding 
coordinate  system.  Figure  2  exhibits  the  definitions  of  the  angular 
measures  of  the  rotor. 


The  kinematic  boundary  conditions  on  all  interacting  lifting  surfaces 
expressing  the  impermeability  of  the  boundaries  can  be  written  in  the  gen¬ 
eral  form  as 


©  ©  © 

WS  *  II  LRKRSdSR  +  II  LDKDSdSD  +  II  LSKSSdSS  ^ 

SR  SD  SS 


0 


© 


WD  =  II  LRKRDdSR  +  II  LDKDDdSD  +  H  LSKSDdS 


(3) 


where  subscripts  R,  S,  D,  refer  to  rotor,  stator,  duct  lifting  surfaces, 
respectively. 


The  kernel  function  K.j  represents  the  induced  velocity  on 
element  j  due  to  an  oscillating  load  L.  of  unit  amplitude  on  ele¬ 
ment  i  .  The  kernel  function  K..  is  the  self-induced  velocity  at  a 

JJ 

point  of  the  particular  lifting  surface  due  to  unit  load  at  each  and 
every  point  on  the  same  surface.  The  kernels  with  two  different  sub¬ 
scripts  represent  the  interaction  effects  from  neighboring  surfaces.  The 
integrations  on  surfaces  S^,  S^,  and  Sp,  are  over  the  rotor  blades,  the 
stator  vanes  and  the  enshrouding  nozzle,  respectively. 


The  terms  Wj  on  the  left-hand  (L-H)  side  of  the  equations  are  the 
known  velocity  distributions  normal  to  the  lifting  surfaces,  nondimension- 
alized  by  the  free  stream  velocity  U.  The  velocities  normal  to  the 
respective  lifting  surfaces  are  the  perturbations  from  the  basic  flow  due 
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to  nonuniformity  of  the  flow  field  (wake),  camber,  incident  flow,  and 
thicknesses  of  the  respective  lifting  surfaces.  In  the  linear  theory, 
their  effects  can  simply  be  added. 

We  consider  two  basic  flows:  a)  one  generated  from  the  hull  wake  and 
measured  in  the  plane  of  the  stator  in  the  absence  of  all  interacting  sur¬ 
faces,  and  b)  the  other  generated  by  the  presence  of  the  hull  and  stator 
together,  measured  at  the  plane  of  the  rotor  in  the  absence  of  duct  and 
rotor.  Thus,  any  harmonic  content  of  the  viscous  and  potential  wake 
generated  by  the  presence  of  the  hull  and  the  stator  will  be  included  as 
an  input  to  the  interaction  problem.  (See  Note  at  end  of  this  section.) 

The  flow  disturbances  considered  in  the  present  study  are; 

1)  The  basic  flows  (hull  wake  and  combinations  of  the  hull  and 
stator  wakes)  both  of  which  will  affect  the  steady  and  unsteady 
loadings  of  all  interacting  lifting  surfaces.  In  fact,  the 
former  will  be  utilized  to  calculate  the  steady  and  unsteady 
loadings  on  the  stator  and  the  latter  will  be  used  to  determine 
the  loadings  on  the  rotor  and  enshrouding  nozzle,  as  will  be 
demonstrated  later  on  in  the  development. 

2)  The  thickness  distributions  of  all  lifting  surfaces  affect,  in 
principle,  both  steady  and  unsteady  loadings  of  the  interacting 
surfaces  as  will  be  seen  in  the  analysis.  These  effects  some¬ 
times  are  omitted  because  of  the  presence  of  the  axi symmet ri c 
duct  conf i gurat ion  and  sometimes  because  the  effect  is  very 
small  in  magnitude,  e.g.,  being  at  the  blade-blade  crossing  fre¬ 
quency. 

3)  The  camber  and  flow  angle  (i.e.,  incident  angle)  of  the  respec¬ 
tive  surfaces  will  affect  their  steady-state  loadings  only. 

Thus,  WD,  the  velocity  normal  to  the  rotor,  is  due  to  basic  flow  dis- 
turbances  in  the  presence  of  hull  and  stator  wakes,  which  affect  both 
steady  and  unsteady  loadings;  the  rotor  blade  camber  and  incidence  angle 
affects  only  the  steady  state  rotor  loading  whereas  the  effects  of  duct  and 
stator  thickness  distributions  may  be  present  in  both  steady  and  unsteady 
state  rotor  loadings. 

The  flow  disturbances  W.  are  made  up  of  the  normal  velocities  on 
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the  stator  due  to  the  hull  wake,  stator  blade  camber  and  incidence  angle, 
and  duct  and  rotor  blade  thickness  distributions.  Details  of  these  con¬ 
tributions  will  be  seen  later  on  in  the  development. 

In  the  linearized  version  of  the  interaction  problem,  the  duct  is 
assumed  to  be  a  cylinder  with  zero  coni city  angle  (i.e.,  a  =  0) .  The  flow 
disturbances  are  those  due  to  non-zero  a  (conic  form)  and  to  duct 

camber,  both  of  which  affect  steady-state  duct  loading  only, and  those  due 
to  rotor  and  stator  blade  thicknesses. 


The  surface  integrals  of  Equations  l,  2,  and  3,  are  reduced  to  line 
integrals  by  approximating  the  chordwise  loadings  on  stator,  rotor  and 
duct  by  appropriate  mode  shapes,  as  in  References  1,  5,  and  6.  The  blades 
of  stator  and  rotor  are  divided  into  small  spanwise  strips  and  the  span- 
wise  loading  coefficients  of  the  chordwise  modes  are  assumed  constant  over 
each  small  strip  so  that  only  the  kernels  need  be  integrated  over  the  span. 
The  collocation  method  is  used  together  with  the  generalized  lift  operator 
technique,  as  in  the  references  cited,  to  determine  the  spanwise  loading 
coeff icients.  In  the  case  of  loading  on  the  duct  of  circular  section,  the 
peripheral  loading  is  expressed  in  terms  of  a  Fourier  series  so  that  the 
peripheral  integration  is  easily  performed,  and  the  chordwise  loading 
coefficients  are  obtained  by  the  collocation  and  generalized  lift  operator 
methods . 


The  kernel  functions  are  derived  by  means  of  the  acceleration  po¬ 
tential,  Kdd  as  in  References  1  and  2  for  the  propeller  alone,  and  KftD, 

RR  DK 

L.,  and  K_  ,  as  in  Reference  5  for  the  propeller-duct  interaction.  The 
RD  DD 

kernels  KCD  and  KDC  representing  the  interaction  of  stator  and  rotor  will 
5  K  Ro 

be  developed  following  the  approach  of  References  3  and  k  for  the  counter¬ 
rotating  propeller  system.  The  remaining  kernels  K^,  K^, 
be  derived  following  References  5,  6  and  l. 


.,  and  Kss>  wi 1 1 


The  three  integral  equations  are  solved  by  an  iteration  procedure. 

It  will  be  assumed  at  first  that  duct  and  rotor  have  no  effect  on  stator 

loading  which  will  be  obtained  from  Equation  (2)  by  ignoring  the  first  and 

third  integrals.  On  substituting  that  value  of  in  Equations  (1)  and  (3), 

those  equations  will  be  solved  by  the  iteration  procedure  outlined  in 

References  5  and  6,  thus  obtaining  values  of  LD  and  L_.  The  values  ob- 

tained  for  L  and  L-  are  then  substituted  in  Equation  (2),  which  is  solved 
R  w 


k 


for  a  new  Lg.  The  new  Lg  is  next  used  in  Equations  (1)  and  (3)  which  are 
put  through  the  iteration  process  again.  The  procedure  is  repeated  until 
stabilized  values  are  secured. 

The  f i rst  set  of  i terations  will  yield  first  approximations  of  the 
loadings  by  solving 

ws  ■  II  LsoKss‘Ks 
ss 

WR  '  II  LS0KSRdSS  *  II 

ss  SR 

WD  ”  H  LS0KSDdSS  “  II  LR0KRDdSR  +  II  LD0KDDdSD 

SS  SR  SD 

Second  approximations  of  the  loadings  will  be  obtained  from 

WS  "  II  LR0KRSdSR  ‘  II  LD0KDSdSD  =  H  LSlKssdSS 

SR  SD  SS 

WR  '  IJ  LSlKSRdSS  =  II  LRlKRRdSR  +  II  LDlKDRdSD 

SS  SR  SD 

WD  -  II  LSlKSDdSS  =  II  LRlKRDdSR  +  II  LDlKDDdSD 

SS  SR  SD 

and  so  forth. 


NOTE ;  If  measurements  are  not  available  of  the  flow  generated  by  the 
presence  of  both  hull  and  stator  at  the  plane  of  the  rotor,  in 
the  absence  of  duct  and  rotor,  corrections  to  the  velocity  on 
the  L-H  side  of  Eq.(l)  must  be  introduced  to  take  into  account 
the  effects  on  the  rotor,  which  operates  in  the  race  of  the 
stator,  due  to  both  viscous  8nd  potential  wake  of  the  stator. 
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THE  VELOCITY  DISTRIBUTIONS 
1.  WD,  Normal  to  the  Rotor 

At  q=0,  the  steady-state  velocity  distribution  normal  to  the  R-H  rotor 
on  the  L-H  side  of  Equation  (1),  after  the  lift  operator  of  order  m  has 
been  applied  to  both  sides  of  the  equation,  is  made  up  of 


w(°.«>/r  )  =  W(0’S){r  )  +  W  (0,S) 
R  {  R'  W  '  R'  W° 


R  (r  )  +  W  <°'ffi>(r  ) 
c+f  ^  R'  WDtR  1  R1 


The  wake  component  (nondimensional ized  by  U)  is  derived  from 


(4) 


1,5 


.  (qR,m) 


(qj 


where 


"W 


,<V 


*  s  J  -V  "V 


-,<*RtpB0 


d9 


Of 


(5) 


9, 


RO 


bR 


=  q  -harmonic  of  wake  velocity  normal  to  the  rotor  blade  in 
K 

the  presence  of  the  hull  and  stator 

=  0^-0^  coscp^  ,  angular  position  of  control  point  with  re¬ 
spect  to  midchord-line,  radians 

=  angular  position  of  midchord-line  of  the  projected  blade 
from  the  reference  line  through  the  hub 

=  projected  semichord-length  of  the  blade  in  radians 


§(m)  =  lift  operator  function 


With  I  <m)  (x)  =  £  J  *(m)e  "dcp^  (see  Appendix  A),  the  wake  harmonic  com- 


.  tt  ixcosqs 

1  r  */-•, .  a 


ponent  is  defined  as 


(qB) 


and 


.(qR,m)  Vw  "  “iqR0R 

W  R  fr  1  =  — -  (r  le  K  R  8  } 

W  R;  U  '  R;  '"r  bR' 


v(0) 

gw°’m)(rR)  =  -V  {rR]  |(m)(0) 


(6) 


(7) 


R  i  ght-handed 


jL 
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The  nondimensional  normal  velocity  component  W.  ,  due  to 

Kc+f 

effects  of  camber  and  incident  flow  angle,  which  is  present  only  in  the 

steady  state  (q  =0)  since  the  blades  are  considered  rigid,  is  given  as 
R 

the  sum  WD  +  W  ,  where 
Kc  Kf 


u  (°‘m)  %  ; - 

Rf  (r  )  =  -  +  a 


2  r  2  ; 
R  R 


(8) 


where 


aR  'nrR0/U 

Q  =  magnitude  of  angular  velocity  of  rotor 

r  „  =  radius  of  rotor 
R0 

0  =  geometric  pitch  angle  of  rotor  blade 

PR 

P  =  tan  ^(1/a  r  )  =  hydrodynamic  pitch  angle  of  assumed  helicoidal 
surface 


u  (°»S) 


<rR> 


"CR<rR> 


"  -  °f(rR»SR) 

j  §(m)  - dcp, 

O  R 


a 


(9) 


where 

f(rD*SD)  =  camberline  ordinates  from  the  face  pitch-line 

R  R 

SR  =  ( 1 -cos^) /2  ,  chordwise  location  as  fraction  of  chord 
length  CR 

CR  =  chord  length 


(This  component  is  derived  in  Reference  8  for  arbitrary  camber  shape.) 


The  nondimensional  normal  velocity  component  due  to  the  effect  of 
duct  thickness  on  the  rotor  is  derived  in  Reference  6  for  a  modified 
lenticular  chordwise  section  (see  Figure  3  represented  by 

f(6a)  *  ii  [to-do]sin20crt-  do(l-cos0a)i  ,  0  £  ea  S  n) 


as 
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2R„r„a, 


u-  <0,»)(rR)  ,  .  -0-5-  J  [k  -o0crR)Ko(kRD)] 

1  tt  v  ]+a2  r2  ° 

R  R 

rr  -i  'a  €D^  ,k0. 

•R'P*lL(2to‘2do)F(k)+idoG(k)jc  |(m)(_bRj|dk 


(10) 


where 

Rp  *  radius  of  cylindrical  duct 

t  =  maximum  duct  thickness 

o 

dQ  =  semi -thickness  of  duct  at  trailing  edge 

=*  axial  distance  between  rotor  plane  and  duct  midchord 
(ep  is  positive) 

and  F(k)=[sin(kC0)-kC0cos(kC0)]/(kCD)2 
G(k)  -  sin(kCD)/kCD 
Cp  =  semichord  of  cylindrical  duct 

lQ(  )  and  K  (  )  are  modified  Bessel  functions  (see  Eq .(19)  of  Ref. 6) 


When  qR  j*  0  (unsteady  cases),  the  velocity  distribution  normal  to 
the  rotor  on  the  L-H  of  Eq.(l),  after  the  lift  operator  has  been  applied 
to  both  sides  of  the  equation,  is  made  up  of 


W  R  (rR) 


(11) 


where  the  W w  is  given  by  Eq.(6). 

As  shown  in  Reference  6,  for  an  axisymmetric  duct,  with  d  =  constant 

o 

over  the  circumference,  (as  in  the  pump-jet  system),  there  is  no  effect  of 

duct  thickness  on  the  rotor  when  q  j*  0 . 

R 

As  noted  in  the  preceding  section,  if  the  wake  of  the  stator  has  not 
been  measured,  additional  normal  velocity  components  must  be  included  due 
to  the  potential  and  viscous  effects  on  the  rotor  of  the  race  of  the  stator. 
These  are  derived  in  Appendices  L  and  M  as  suggested  by  Dr.  John  Breslin. 
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2.  Wg,  Normal  to  the  Stator 

In  the  steady  state  (q  =0) ,  the  nondimensional  velocity  distribution 
normal  to  the  stator'on  the  L-H  of  Eq.(2),  after  the  lift  operator  of 
order  m  has  been  applied  to  both  sides  of  the  equation,  is 

y (°»m)  ( r$)  =W^°-m)(rs)  +W$2’m)(rs)  +WD^,m)(r$)  +  WR(°,rn)  ( r$) 


(12) 


Here 


wV0’"^)  =  —  (rs),(m)(0)  (wake  of  hull  alone  in  plane  of  stator) 


(0) 


w'“'S)<rs)-  -Vl+  *|r|  [eps(rs)  -  B  ( r$)]l  <">  (0) 


(12a) 
(1 2b) 


where 


a  = 


-5  rst.„<yrs) 

»  *  6ps<°-7W 


at  r$  -  0.7rR0 


-  (0  ft  ^l+»srs  f  -  af(rS’ss' 

"s'  <rs>  ■  5-c^T  J  *<">  >*s  <% 


(12c) 


(cf.  Eq.(9)  for  details.) 


The  velocity  due  to  the  effect  of  duct  thickness  on  the  stator  can 
be  shown  (see  Ref.  6)  to  be  equal  to 


(0,m) 


DtS 


2Vsas 


TT  V 


>+asrs  ° 


I  lkIo(krS)Ko(kRD>l 


'R'P-{[(2to‘2do)F(k)+idoG(k)]e 


•ikps.  - 

Vas 


e  +e 
D  S 


|(m)(“ri)}dk  O3) 


which  is  Eq.(10)  with  stator  geometry  substituted  for  rotor  geometry.  Note 
the  factor  exp(-ikCg)  which  is  the  result  of  the  substitution  xs“cf>so^aS+€S’ 
where  is  the  axial  distance  of  the  stator  from  the  rotor  and 

*The  stator  has  the  geometry  of  a  left-handed  propeller. 
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cpso  =  a5“®|35COSCPa  (see  Eq . (5) )  -  In  this  case  the  stator  being  forward 
of  the  rotor  is  negative. 

The  component  WR  due  to  the  effect  of  rotor  blade  thickness  on 
the  stator  can  be  showS  (see  Appendix  B)  to  be  given  in  the  steady-state 
by 

I  -  fcu.^ok). 


tt2  Vl+a2  r2  PR  b* 


/°S  aR 

■il— - -  +  e  )u 

Vas  8«  sJ 


I  (ueb  /e  )|dudp 


where 


(IK)  .  r'o!'jp(i)Kot,jrs)  for  pR<rs 

l0(xrs)K0(dPR)  for  r$  <  PR 

F(u,pr)  -  {Sin(uebR/aR)  -  (u9bR/aR)cos (u8bR/aR)}  /u! 


Nd  =  number  of  blades  of  rotor 

In  the  axisymmetric  duct  case,  which  is  presently  under  consideration, 
there  is  no  effect  of  the  duct  thickness  on  the  stator  or  rotor  for  the  un¬ 
steady  flow  case,  i.e.,  q^O  (see  Eq.(2)  of  Ref.  6),  so  that  the  velocity 
distribution  normal  to  the  stator  is 


.  (q«.»n>)  (qc,m)  (XN  ,m) 

WS  (rS)=WW  (rs)+i?RS  (rS)  ’  je  •=  1 ,2,3,  . . .  (15) 


where 


_  (qs,m) 

Wy  (cf.Eq.(6))  is  due  to  the  wake  of  the  hull  only  measured 

at  the  plane  of  the  stator  (15a) 


and  the  effect  of  the  rotor  thickness  (see  Appendix  B)  is 
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_  (IN  ,m)  4a=N  r 

wR  c  R  (O - R  R-S  -- 

1  rr  2  J 1  +az  r2 

S  S 


Vs] 


P»  X 


J  e^“  ~T  (PR)  ^]+  3rPr  J  F(u,pR)[G2(u)-Gs(-u)]dudpR 
u "  o 

(X  =  +l,+2,  ...  )  (15b) 


PR  bR 


GR(u)-'iN1!(lU-,aR%ll>R)KiNB(lU+‘R"Rlrs)LaSU+%(vR  '  T')] 


’R . .  ■  R 

iu(°S/aS-°R/aR  +CS> 


■  (">[(-%(-  J)  *  ^Xs] 


3.  WD>  Normal  to  the  Duct 

In  the  steady  state  (q=X=0) ,  the  nond i mens i ona 1 i zed  velocity  dis¬ 
tribution  normal  to  the  duct  on  the  L-H  of  Eq.(3),  after  the  lift  operator 
of  order  m  has  been  applied  to  both  sides  of  the  equation,  is 

w(0,m)  (m),  ,  a  {0,m)  rs  (0,m)  g  (0,m)  (16) 

D  1  W  WDC  RtD  StD 

where  the  first  component  is  due  to  conicity  angle  a,  the  second  to  duct 
camber,  the  third  to  rotor  blade  thickness,  and  the  fourth  to  stator 
blade  thickness. 

It  should  be  noted  that  in  the  linearized  version  of  the  interaction 
problem,  with  the  duct  assumed  to  be  a  cylinder  with  no  conicity  angle, 
it  is  assumed  that  there  is  no  contribution  to  the  normal  velocity  on  the 
duct  surface  (i.e.,  in  the  radial  direction)  due  to  the  wake  generated  in 
the  presence  of  both  hull  and  stator  together. 

Reference  6  shows  that  for  axisymmetric  ducts,  and  assuming  a  modi¬ 
fied  lenticular  camber  distribution,  namely, 

«(?„)  *  (mx  +  ^)sin2<pa  "  I  <V,"coscp!P  *  0  s  <Pa  S  n 
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where  m  is  maximum  camber  and  d  the  semi -thi ckness  of  the  trailing 
x  o  J 

edge  of  the  duct,  the  component  due  to  duct  camber  (see  Eq.5^  of  Ref. 6) 
is : 


aD°,m>  =  r  {(2mx +  d0} 1 1 (m) <°>  -  t  1 (m)  (0)} 


07) 


(l/m>(0)  is  defined  in  Appendix  A.) 

The  same  reference  derives  the  velocity  component  due  to  the  effect 
of  rotor  (propeller)  thickness  on  the  duct  (see  Eq . ^3  of  Ref. 6)  as 

^a«N. 


;  (o,S)  .  _  f  J^ST:  >.  1°  (p  , 

R  D  __2  w  $  a.  t;  C  vrR7 

t  "  '  url  '  ‘ 


00  ,  -  iu  (e.-C  )  aR  , 

J  ulo(upR)K,,(^RD;r  '  ,0,pR).|m  Partje  Vm\uCD)]dudpf 


(18) 


where 


F(u,0,Pr) 


sin(uebR/aR)  -  (u0ho/ap)cos(ueho/ao) 


bR  R' 


bR  R'' 


By  analogy  with  the  above,  the  component  due  to  the  effect  of  stator 
thickness  on  the  duct  would  be 


4a*  N 


"s'S’i)  -  -  f  ^  t  bj 


"  r 

J  ulo(ups)K, (uRd)f(u,0,p$) *lm  Partji 


“iu(CD~Cs"aS/aS) « 


^uCpjjdudp, 


(19) 


where 


F(u,0,ps)  = 


_  sin(u9bS/aS)  ‘  ^uebS/aS)c0s(u6bS/aS) 


When  ,  unsteady  flow,  Reference  6  gives  the  nond i mens ional i zed 

<XNR,m)  i£NRcp  i ^N-Qt 

velocity  component  as  WD  .  e  e  *  ,  where 

RtD 
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(4NR,m)  iZ*aoNR  - - -  pD  t 

RtD  J  V,+aRpR  e - r  (pr) 


•|uliNB<uVK>N|i(“V<!  "1* 


R  faR 

Ju(VaR/aR) 


l(m)(-uCD)F(u,-XNR,pR) 


_e-|u(VWl(m) 


and 


F( 


•  m  (uCD)F(u,£NR,pR)]dudpR  (20) 

U./N..P  )-  sin( 6bR/aR)-(“-aR^ J(6ba/^)cos( (u-aBmB) ehg/.J 

M5  ~ 


A  similar  formula  can  be  derived  for  which  however  can  only  be 

effective  when  i  =NR  (blade  crossing  frequency),  since  N$  is  usually 
not  an  integer  multiple  of  NR  ,  and  thus  is  negligibly  small. 
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COMPONENTS  OF  THE  SYSTEM  OF  INTEGRAL  EQUATIONS 


1)  Kernel  Function  K 


RR 


From  References  1  and  2,  the  first  integral  of  Eq.(l)  can  be  shown 

to  be  equivalent  for  each  frequency  qD  to 

K 


I,  =  e 


I  S5|  n5,  Lr  (PR)KRR  (rR’CpR0’PR’eR„:qR)dpF 


PR 


where 


-iq  Aa 

_  (m, n) 
Knn 

“nr 

I -  2 

RR 

'•’"V^ro 

•aR^,+aRrR  "l^OrAnr 

9, (u) -g,f-u) 


(21) 


where 


g.(u)  =  ( I K)  B (u) e 


•  U  A 

i  —  La 
aR 


("<)„,  ■  { 


VIU+VlVPR>Vl',+Vl,VrR)  for  pR  <  rR 


lm1<lu+aRilNFtlrR>Km1(lU+an'lNRl,,F!)  f°r  rR  PR 


"1 


B(u) 


(aRU+aR£lNR  +  T5  ]  (aRU+3R£lNR  +  7  ) 


•  a^ebr,A(">((V  a^6b  ' 


p^.  =  fluid  mass  density,  slugs/ft3 


rR0  = 
La  =  ar-ap  * 


n  = 


rotor  radius,  ft 

difference  between  skewness  of  the  blade  at  control 
point  r  and  skewness  at  a  loading  point  p  ,  radians 

Q  i-  /U  and  p  and  r  are  also  nondimensional ized  by  r 
RO  o 

angular  velocity  of  rotor,  radians/sec,  U“  freestream 

velocity,  ft/sec 
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r  p 

0^,0^  -  subtended  angle  of  projected  semichord  of  blade  at  r, 
at  p,  radians 

lm(  ),«„,(  )  modified  Bessel  functions  of  the  first  and  second  kind 
and  **  0,  ±1 ,  ±2,  ... 

(qR*n) 

The  L0  (pD)  are  the  unknown  spanwise  loading  coefficients  after 

R  R  (q0) 

approximation  of  the  unknown  loading  function  L  R  (PD»®Dr>)  'n  the 

R  RO 

chordwise  direction  by  Birnbaum  mode  shape 

(qp)  *  (qR.n) 

L  V  -  -S,  L  ' 

n=l 

ir  (^.0  ea  •  (qR»n)  _  , 

=  w{LR  (PR)cotT+EL  (pR)sin(n-l)eJ 

n=2 

where  0Rg  =  op-e^cos0ft)  and  after  the  subsequent  chordwise  integration 


over  0 


RO 


JbW5Ua 


■  i  -iycosB-, 

My)  “  ^  J  ®(n)e  sin0adea 

o 


(See  Appendix  A.) 

The  cp  (=ar-0^"coscp  )  dependence  is  eliminated  by  operating  on  both 
RU  D  O' 

sides  of  the  integral  equation  by  the  “generalized"  lift  operators  $(m). 
The  factor  I (x)  in  the  kernel  function  is  the  result  of  this: 

i«w 

o 

(See  Appendix  A.)  Equation  (21)  has  an  integrable  singularity  at  u»0,  the 
value  of  which  is  determined  by  (.’Hospital's  rule  as  shown  in  Appendix  C. 


2)  Kernel  Function  K 


DR 


When  the  control  point  is  on  the  rotor  and  the  loading  point  is  on 
the  cylindrical  duct,  the  induced  velocity,  nond i mens i ona 1 i zed  by  U  , 
normal  to  the  rotor  blades,  is 

'2"  p  Q  JI  LD^^D’PD’Ve  *  KDR(XR’rR,9R0;?D,PD,6D: VdSD 


or 
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v  Zr  r  (*a)  !^nK  „  ^  j p 

l2  =  I  ,1  LD  ®  KDRPDd8Dd5D 

\2=°  o  2C0 

where  ^  ^  -p^  =  duct  loading,  lb/ft  (see  Reference  5) 


(22) 


KDR~ 


limit 


xR^R0/aR  »  0  - 

pd"*rd 


„  „  Xr'?D  iX2aR(T*XR+V 

o=re 

sr  sr  j  — - - dT 


DR 


( 


d 

‘*S 


-  r~) 

rR  **> 


3n„  Op 

D  D 


"dr  ■  +  rR  +  p0  •  ZrRpD  “s(teD  -  ?R0  +  “>}" 


CD„„  =  CTn  "  ®unCOSW  ,  0  ^  CD  S  TT 

RO  R  bR  ~a  ra 


The  loading  will  be  expressed  in  a  Fourier  series  as 


LnM(5D,PD,eD)  =  2  L^’^(U. 


•ili0r 


(23) 


at  Pq=Rq*  The  reciprocal  Descartes  distance  1/R  can  be  expanded  in  the 
form 


im, (+9_-cp  +nt)  ob  . 

„  2  °  «  I  (24) 

DR  m2=-®  -os 


1  I  f 

- —  =  —  Li  e 


(rR  <  pD  in  the  limit  as  pQ  -  f%  .) 


From  the  ©^-integration,  it  is  determined  that  m2  =  4  ,  since 

{  (25) 


2n  i(m2-n)80  ^  j  2n  for  m2  -  =  0 


I  e 

o 


d9 


^  0  otherwise 
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Also,  since  the  L-H  of  Eq.(l)  is  an  exp(iqRfit)  function  of  time  and  l2 
is  an  exp[ i (\2+m2 )0t]  function, 


^2  +  m2  =  qR 


(26) 


where  q„  s  0,  ^  0.  The  double  series  can  thus  be  reduced  to  a  single 

infinite  series.  Equation  (22)  becomes 

_  j)  2e 

1  2 

A.  2= 

R  "*  rD  "D 


limit 


GO 

1  —  S 

x2=0  4ttp  u2  xR-*CpRc/aR  2c 
ni2=qD-\2  p  ”* R 


»  (A2,m2)  8  8 

J  ld  <  V  3^  3^ 


-'Vno  ''VWVV  Y'° 

•e  e  J 


-'WV  J  '  (IMrR)K  (Iklp,,)  eiTWd5[ 

-co  z  Z 


(27) 


After  the  r- integration  is  performed  and  the  derivatives  and  limits 

C  T 

are  taken,  the  generalized  lift  operators'  are  applied.  Equation  (2) 
becomes 

"o  J  «Sr,a>-«D 

m  I  a.2  u  Z  C  q 

where  K^R2,m^  is  the  modified  kernel  after  the  cpQ- integration: 


:(rm  ,m)  _  _ I 


■,rn2aR 


DR 


R  R 


m2- 


iTlaR,m2'qR,LaR(m2-V+7j,m2(aRlVqRlrR^LS-l(aRU2'qR,RD) 


+Km,+1 ^R^'V  V  J 


1  -laR(m2-qRHWaR)  .  (m) 


'W(qRV 


x  -Jk(§n-  ^) 


+a. 


k|kllmJlk|rR)[l<ms-l(lklRD)+Km2-H(|klRD)]|(m)((m2-fe)ebR)e  °  ^  dk 


k  -  aR(m2  -  qR) 


(Cont 'd) 
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m_  » 

-Jj 

'  p»  -  CD 


k  -  aR(»2  -  1r) 


(28) 


evaluated  at  m2=clR"^2  *  The  rotor  radius  rR0  is  !ntroduced  in  the 

denominator  of  the  f i ret  factor  because  now  rR,  R^,  and  are  fractions 

of  r  and  a  =  r  SI/ U  .  The  k-integrals  have  integrable  Cauchy-type 
RO  R  RQ 

singularities.  There  is  no  other  singularity,  as  can  be  seen  from  the 
original  kernel  of  Zq.(22)  since  rR  is  always  less  than  P^Rq. 

If  the  chordwise  loading  on  the  duct  is  approximated  by  the  Birnbaum 
mode  shapes 


0  00  — 

L(X2,m2)(  )  _  1  ^(Xs.rte.l)  cot  ^  +  _S  A(Xs,m2,n)sin(n-l)ea  }  (29) 


n=2 


where  6a  is  defined  by  ?D  «  eD  -  CQ  cosD^  ,  0  s  0ff  s  n  (see  Figure  1),  then 
the  integration  over  §D  is  easily  accomplished. 


a  Tt 


(X2>m2,n)Qf-^-  (m2,m)j,  sjnQ 


I  L^’rn2)(VRDR2,m)dV  W  -S  I  AU2>m2’n)®(n)KDJn 

2CD  n=1  o 

V  *(X2,m2,n)-  (m2tm,n) 

A  ^DR 

n=l 

where  0(n)  are  the  chordwise  mode  shapes  given  in  Eq.(29), 
and  A(Xs’rn2,")  =  cDA(X2’m2’n) 


a  or 


(30) 


and  ,rn>  n) - 1 


■im2oR 


^Pf^RO  * +  *1  i 


lTTaR|m8'qRl[aRK-qR)+ 7  ] 

rR 

1  w  ’  ~R  R 

|m2(aR|m2-qRirR)[Kms_,(aR|m2-qR|RD)+Km2+,(aR|ma-qR|RD)] 

-iaR(m2-qR)(eD-oR/aR)  (-}  (-} 

•  e  1  'qR6bR  A  — aR  qR  °D 


-^(•oST5 


1kl'm  (lklrR)[Km  -lOklRj^  Rdk 


rR 

Then  letting  u  =  k  -  aR(m2-qR) 


-(m2=qR-X2,m,n)  i 


'DR 


k  -  aR(m2  -  qR) 


rR  /^r/’V^R^V 


(31) 


4npfu  rR0 
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•  {ing2(0)  +  J[g2(u)  -g2(-u)]^} 


(31a) 


where 


r  2  "Un  r  ) 

92<u>-LV+aR<'W  +  4\  Llu  4  V»2  -  V'J  '  * 

*  n 


lmsOu+aR(,V‘:|R>lrR)  [Km2-lOuaVVVlRD)a*m2+l(lu<VVVlRl>)] 


See  Appendix  D  for  the  evaluation  of  the  singular  part  of  Eq . (3 la) . 


3)  Kernel  Function  K 


When  the  control  point  is  at  (xR.rR»tPR)  on  the  rotor  and  the  loading 
points  are  at  (Sg»Ps*®gg)  on  the  stator,  the  nondimensional  induced  velocity 
normal  to  the  rotor  is, 2,3 


"  r  r  (Aal  U,Ot  , 

n  ls  «vva  s  «s 


^3  -0  o  Pc 


cd  ~  Z,  °  sinB  d0  dp_ 
SR  a$ps  a  a  KS 


where 


eso  ■  °s  •  ebs“se« 

_  2S  I'm  »  [■“  l'l3LaR(T,'V"5Sn]  i  {  ]  \ 

5R  ^pfUa  nt,  **  L  6 


rsrm  {(t'  •  !s)2  4  **R  *  "l  ■  2Vscos(9S  •  0}' 


■  {(’'  -  HJ 


S'  *  rR  *  PS  -  2rRPSCO5(~9S0-tPR0+  nt 


-  U‘ 


and  6$R  -  0  means  that  xR  -  cpRQ/aR  and  g£  -  6S0/as  +  cs 
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The  inverse  Descartes  distance  1 /R  is  expanded  as 

5  R 


1.1?-  lm3SSR  *.  ''<T'-Wk 


SR  m3— » 


s  '»3<lkl‘,s>lWlkl'Vdk  <33> 


(for  „s  <  rR)  with  PSR  =  -6S0  -  <fR0  +  fit  -  6^ 


The  derivative 


"s 


(  a  i  a  \ 

^rq-7%^-0) 


is  the  directional  derivative  normal  to  the  stator  blade. 

After  the  r'-integration  and  derivatives  are  taken,  the  kernel 

function  K  becomes 
5  R 


K  ,  K  - *-*- 


SR  SR  asp$ 


RR 


Ms 

00 

E 

E  e 

n=l  i 

m3  =  -ee 

m. 

,  -ia„\. 

3 

L  R  . 

_kR 

r 

.  .  (9ro  eso  \ 


'm3(aRX3ps)|WaRX3,V 

'Port  Gc 


.  */  m,\  /  m, \  m, 

-iv-  4)(v  4) 


1  (Iklp  )K  (Iklr  )  ik(— 52.  -  _§0._c  ) 

m,\  rru  S  R  »  3  _  S' 


"3 _ "  '”R  as 


k  +  aR\3 


dk} 


From  the  time-dependent  factors  on  both  sides  of  Eq.(l) 
X3  +  m3  “  qR 

and  from  the  summation  over  the  stator  blades 

Nc  •/,  N.  for  4  N 

S  -i (\,+m,)6_  r  S  3  3  3  S 

s  e  3  3 Sn  -  \ 

n=l  0  otherwise 


(34) 


X3+W  Vs  '  l3  *  0 
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Thus  the  frequencies  of  the  first  equation  are  limited  to  zero  and 

positive  multiples  of  the  number  of  stator  blades. 

Since  =  q^  -m^  ^  0,  £  q^  and  the  double  series  over  X^  and  m^ 

can  be  reduced  to  a  singl  e  infinite  series. 

The  unknown  loading  function  is  approximated  as  before  in 

chordwise  direction  by  Birnbaum  mode  shapes.  After  the  chordwise  integra¬ 
tion  over  ©a  and  application  of  the  generalized  lift  operator,  can  be 
wri  tten 


s  i  ££  L'x-"> 

5  \j=0  Ps  m=l  n=l  S 


3’  (PS)RSRm,n)(m3=VX3)dpS 


where  the  modified  kernel  is 


RSr’  ^  ^m3~  qR‘V  I- 


^nPf U2  rR0  as  V,+a2r2 

K  R 


{LaSSR<qR-“3)+  i  ]  :!]  V^rV^ ’sWVV’V’J 


rR  3 


-im  a  “'Vo  iaR(qR'ml)(W/a5)  fnU  aR  \  (m\ 

•e  3Se  3  A(  )((^(Vm3)-m3)ebS),()(qRV 

j  *7  m3V  m3'\  lm3(*klPs)Km3(lk|r  )  -ike  ("V  a^)CTR  ~l(m3+^aS 

•rf  JL\ask'  4Xv+  $)  L  *  e 

PS  rR  k  +  aR^R-m3) 


where  now,  a,k,r,  and  p  are  nondimensional ized  with  respect  to  rD_. 

RO 

Let  u=k+a^(q^-mj) .  The  kernel  may  be  written  as 


SR'")(m3'VX3)'{- 


^PfU2r  0  ac  'vM+aV 
t  RO  s  R  R 


/'Vr  eiaR<V"3>(V^>J 


.  •  g,(u)  -  S3(_u)  \ 

•  3  u 
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where 


m3“i"  m,3  -j 

93  u  °  LVsVW-  7iV-aR<V"3>  *  f  J 


CR  °S 

“iu(eS_  7"  +  a  5 

•|m3(lu-aR(^R-m3)lps)Km3(lu-aR(qR-m3)lrR)  e  R  b 

for  Pg  <  rR.  Note  that  for  rR  <  p  ,  these  are  interchanged  in  the  modified 
Bessel  functions. 

See  Appendix  E  for  the  evaluation  of  the  singularity  of  KCD  as  u  -»  0. 

5  R 


4)  Kernel  Function  K 


RS 


If  the  control  point  is  at  (x',r  ,cp .)  on  the  stator  and  the  loading 

points  are  at  (§R,  pR,0R)  on  the  rotor,  then  following  the  development 

for  the  other  kernel  functions  it  can  be  shown  that  the  nondimensional 

induced  velocity  at  the  control  point  due  to  N.-blades  of  the  rotor  will  be 

R 

given  by 


os  TT 


I.-  E  r  f  L a*}(p  ,e  )e'^QtKl  ^Sine  d6  dpR 
4  .  _n  J  J  lR  vVVe  *RS  a  «  a  ofR 

V°  o  pR 


RS  Vr 


(38) 


whe  re 


krs= 

,  nr 

1  v>  1  i  m 

4TTPfu2  n=l  6RS"*° 

8  ?  iXl,taR 

5nJ  J  e 

S  -00 

a 

=  -rS  (a  l 

.  J-l_' ) 

°"s  ' 

^  i  +as  rs  * 

1  V 

X  1  = 

xs 

*so/as  +  €s  (es 

negative) 

1  nt 

r^  =  0.7  radius 

as  = 

rStan0PS(rS) 

t —  l  TT — j  dT 1 


RRS-{(T'-VS+  1+  [’R-2rsVOa[6R04<PS0-nt+5Rn-aR!''  '  *s>  ^ 


and  by  6R$- 0  is  meant  x'-cp^  +  eg  and  CR  -  6R()/aR  . 
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Expanding  the  Descartes  distance, 

1  -I  ?  „'m«eRS  7  '<T'-Vk 


J _ =1  E  e 

R  TT  "  e 

RS  m4=-» 


I 


'n,4(|k,^)S(,k,rS)dk 


(39) 


for  p  <  r  ,  otherwise  pp  and  r  are  interchanged.  Here 
Kb  b 

eRS  =  0RO  +  ^SO  -  Qt  +  eRn  -  aR(T'-x^)  . 

After  performing  the  T‘-integration  and  taking  the  derivatives  in  the 
proper  order,  and  taking  the  limit 


+anPo 

R  R  _ 


KRS  KRS  “  aRpR 


NR  -  ■ 

E  E  e 


■im4nt  i(m4-X4)eRn  in,4(eR0-KpS0) 


4npfU2aR  Vi+a|r2  n-’ 
b  b 


9, 


iaRK‘V  (if  +  es~  if) 


9, 


R0\ 


mi. 


*  ILVrK"^)"  7'jLaR(m4"X4)+  7J 


41. 


•  ,m4(aRlm4  '  X4j  WaR,m4  '  ^,rs) 

Vcn  ®Rn 

l*'  "W  m4f  m4(lk,PR)Km4(,k,rS)  ,k(T+  V  aR  >J, 

-  tt  Hask-  ~AaRk+  - - - -  e  dk 


k-aR(m4-X4) 


(40) 


The  time-dependent  factor  on  the  L-H  of  Eq.(2)  is  exp(iq^nRt)  and  the 
time-dependent  factor  of  l4  on  the  R-H  side  is  expL i (\4~m4)Qt] ,  therefore 
q$  =  X4-m4> 

NR  i(m4-X4)6Rn  ,  Nr  for  (m4~\4)  =  \NR,  \  =  0,  ±1 ,  ... 

Also  2  e 

n=l  0  otherwise 


so  that 

q$  B  i4NR,  \  -  0,  +1,  +2, 


and  X4-m4  *  0  . 


23 


After  the  chordwise  integration  over  0Q  is  performed,  by  represent¬ 
ing  the  chordwise  loading  distribution  by  the  appropriate  mode  shapes 
8(n)  (see  ^ )  and  the  generalized  lift  operator  §(m)  is  applied,  the 
integral  l4  becomes  for  each  q^,  m  and  n 


>4  = 


iqsOt 


CD 

£ 

v° 


L 


(^4,n) 

R  <VK 


(m,n) 

RS 


dpR 


(41) 


where  the  modified  kernel  is 

N, 


K 


(m,n) 


RS 


(m4  =  VqS}  =  ' 


im4^c7R+°S) 


^Vro 


V  4V  2  4\ 

•ivasaRV7fAaRqs-7> 


mIA  "iaRqS(eS+i7*^) 


,m4^aRqSPR^Kn’4^aRqSrS^ 


•  A<n)(H +  VeJ  (i-VrVeJ 


as  ctr 

m,.\ /  m,.  \  ik(eS+i7’^ 

„  -"S  '  3/V'R"’  p2  /  6  PR)Kmi^(lk' rs) 

rS  R 


-  * !  (v-  SX-t*  £ ) 


^(H  -  I^r)  1  <S)  (<"V  r>6bs)  k+a75~  }  W 


S  *'*"  k+aRqS 


Let  u  =  k  +  a  q  ,  then 
K  5 


_(m,n) 

KRS  (m4  “  VqS> 


rS  ’WV 


4npfUSaRrR0 


s  s 


c  o 

-ia  q  (e.+  -  -£) 

RMSV  S  a.  a  ' 


s  V  W.).lj  Vl^ldu} 


where 
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a,  a 

(  mi4\  /  s  m4\ 

9i4(u)=  ^asu-asaRqs-  — AVW  ^ 

rS  PR 

'  lm1,(lu-VslpR)l<™^-Vslrs) 

•  A(n,(<w  ir>eb»)l<m)((',ni. 

R 

(for  pR  <  r$)  .  See  Appendix  F  for  the  singularity  of  KR$  at  u  *  0. 

5)  Kernel  Function 

When  the  control  point  is  on  the  stator  and  the  loading  point  is  on 
the  cylindrical  duct,  the  nondimens ional  induced  velocity  normal  to  the 
stator  blades,  the  second  integral  of  Eq.(2),  is 


v  ±  ^ 


Ts  %  ■  t>ebs) 


(43) 


'5*  £  U  ldW 

X5=0  o  2cd 


where 


1  limit 


DS 


4ttP  U2  xs^S/aS+eS  unS  “PD  -« 

'  V  R0 


a  a  s  ’D 


R0S  =  {T2  +  1  +  PD  ■  2rsPDcos(eo-,l,S)} 


(44) 


dT 


DS 


■  {^  +  1  +  PD  ’  ^sV^D^V}' 


The  loading  will  be  expressed  in  a  Fourier  series  as 

Us) 


l-D-'tWV  ■  2  L'Xs’1‘>(SD)e",‘eD 

^BbQD 


(45) 


at  Pg  =  R^.  The  reciprocal  of  the  Descartes  distance  can  be  expanded  in 
the  form 


J_.i  £  .'WW  j  ,  ,|k|  K  (tkle  )eiTkdk 

*D$  ms—' »  -®  5  5 

since  r$  s  p^  in  the  limit  as  pp  -  Rp  . 
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From  the  0^-integrat ion  it  is  determined  that  because 


f  ei(m5^)6D  deD  =  { 


2n  for  m5-(i=0 
0  otherwise 


(47) 


Since  the  L-H  side  of  Eq.(2)  is  an  exp(iq^QRt)  function  of  time  and  1^  is 
an  exp(i\sQt)  function  of  time 


b  '  qs  ■  l\ 


m 


With  the  substitution  of  (45 , (*+6) , (47)  and  (48),  and  after  the  T-integration 
and  the  derivatives  and  limits  are  taken  and  the  generalized  lift  operator 
is  applied  using  the  complete  orthogonal  set  of  functions  designated  as 
$(m)  ,  becomes  for  each  m,  order  of  lift  operator, 


-  -  -  (qc.mr)  iq«.nt  _  (mq,m) 

'i'  =  L0  5  Ss5  «» 

qs~0  2Cd 


where  the  modified  kernel  (after  the  cp^- integration)  is 


_  (m5,m)  1 


.'m5as 


DS 


4irp  U2  r  V 1  +&t  r2 

yf  RO  S  S 


(49) 


r  .  m  ,  iqSaR< W  i7> 

l-i"aRqs(.8SaR%  +  7J )  a  '™5(aRqSrS) 


[Km5-l(aRqSR0)  +  Km5+I(aRqSR0)]|(”>((-”5  +^qS)6bs) 


-ik(  W 


lm5<lkl'-s>LKm5-lf|klR»)+,!Vl(lklRD 


(">(K  - 

k+aRqS 


dk 


} 


(50) 


The  expansion  scheme  has  introduced  an  integrable  Cauchy-type  singularity 
in  the  k-integrals.  There  is  no  other  singularity. 
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If  the  chordwise  loading  on  the  duct  is  approximated  by  the  Birnbaum 
mode  shapes  as  in 


(Qc  »mr)  1  (  (qc.'T'c.O  ®, 

s  ■>  (p  )  =  JL  ■//>  *  -> 

■n  '  -n  L 


(qs,m5,n)  _ 


A  w  "  cot^+2  a'  5  5  sin(n-l)6  } 

n-2 


(see  Eq.(29))  then  the  integration  over  is  easily  accomplished. 

„  (qc.mr)  .K.™)  i  “  vr 

J  ko  S( VSs  “V  «  =  j 


n  (q  ,mc , n)  (m-.m) 

*A  S  5  ®(n)KD$5  CDsin6odea 


®  (qc ,mc,n)  (m-.m.n) 

=  2  A  5  *  K  * 


n=l 


DS 


where  A 

and 


E  (*••)  =  CpA^"5 


(see  Eq. (30)) 


(51) 


(m-.m, n)  j 

1/  ?  =  - - - 


'DS 


r$  _  _  e'm5aS 


1 


■  {-’"Vs^sVs  *  -$)  e 


...  'VR<wr> 

S  P 


'  [K™5-l(VsRD)  +  Km5+l<Vs,1D)]|W(<-"5 


-,k<vv 


<lklrS)LV-l(lkl’>0,+K.c+l(lklRD)] 


1  m j  x  ’  ”*  S'  v.  nij- »  u 


-^Xs) 

k+Vs 


dk 


} 


(52) 
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Let  u  =  k  +  aRq$ 


(mr  »  m  »  n)  | 

kds  = - 


rS  iys  iaRqS(Ves"  a^ 

e  e 


^Pfu2rRO  V,+aH 

•  {ing5(0)  +  j  Cg5(u)  -  g5(-u)]  } 


(52a) 


where 


95(u) 


(  m5  \ 

■  ,u  •  aRV  (v  '  Vi"!  s) 


-|u(vv  T1 


Vslrs)LV'(,U  '  Vs^b)  *  V+l(,u  -  VsK)] 


'  1  (m)((-  5^  +  5^  %  -  m5)ebs)  A<n)((-U  +  vOcd) 

Equation  (52)  has  an  integrable  singularity  at  k=  -aRqs.  T^e  value  of 
the  integrand  at  that  point  is  determined  by  means  of  L'Hospital's  rule  as 
shown  in  Appendix  G. 


6)  Kernel  Function 

The  third  integral  of  Eq.(2),  the  nondimensional  self-induced  velocity 
at  a  point  (*5 . rs .9S0)  on  the  stator  due  t0  the  ,oadin9  at  points  (5^, P$ . ©SQ) 
of  all  N$  blades  of  the  stator,  is  given  as 


00  ( x6) 

,6-S  Ls  (Ps.e„). 
xfi=0  Ps 


ix6nt 


*ss 


asps 


sinVeadpS 


(53) 


28 


TR-2173 


where 


Kls=- 


limit  3 


^np f\?  n-t  6SS~° 


U  I 


S  iX6[aR(T'-^)-etJ 


sJ'“snJ  a  /  n 

^  'RSS' 


dT1 


RSS  ~  {<'r'-5S>3  +  rs  +  PS  2rSPSCOS^"6S0  +  ^SO  ■  ®Sr|)} 


6ss~  0  means  that  xs  -  cpso/as  +  es  and  H  ~  eso/as  +  es 


The  inverse  Descartes  distance  is  expanded  as 


i  “  "S  '  (T'“§c)k 

•  J  \u"*>WMr*u 

6 


(54) 


(for  <  r^,  otherwise  and  r^  are  interchanged  in  the  modified  Bessel 

functions)  with  Pec®  ®cn  "  'Pen  +  *  The  summation  over  the  blades  becomes 

bb  bu  bu  bn 


pS  i(m6-X6)eSn  |NS  for  m6‘X6  =  X6NS '  £6  =  °'±1,±2»  *** 


E  e 

n=l 


0  otherwise 


(55) 


Also,  since  the  L-H  side  of  Eq.(2)  is  an  exp(iqsOt)  function  of  time  and 
lg  is  an  exp(iXgfit)  function  of  time, 


Xg  =  q^  =  ,  i^=0,  +  l,+2,  .. 


(56) 


With  these  substitutions,  after  the  T’-integration  and  the  derivatives  are 
taken  in  the  proper  order,  the  kernel  function  =  Kssv,+asps/asps  be‘ 
comes 


Kss" 


“•"WrO  "/Ha|r| 


1  imi  t 

6ss-° 


Im6^0SO-'PSO) 


m6"V*6HS 


[Cont'd] 
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*  KvrV  7)(asaRqs  ‘  7) 


-'qraR(xc“^) 

4)«  VWs’VVs'V 


“  (  ”6V  m6'|  ik<xs"5s'  l">6^k'ps'Km6^k'rs'  1 

"  I  (.ask  +  7'Aask  +  7J e  - j-77- - dk} 

-«  S  PS  k  RqS 


where  a,  k,  r  and  p  are  nondimens ionalized  by  rR^  the  rotor  radius. 

After  taking  the  limit  and  substituting  0^=  ”®bSCOSeo  and  ^SO  = 
-ej  coscp  (0°  skew),  the  chordwi.se  integration  over  cp^  can  be  performed  by 
representing  the  chordwise  loading  distribution  by  the  appropriate  mode 
shapes  0(n)  (see  J^)  and  the  generalized  lift  operator  §(m)  can  be  applied. 
The  integral  becomes  for  each  q^,  m  and  n 

.  (qc.n)  iq^t  _(S,n) 

'6  *  I  4  (ps>‘  KSS  dpS  (5E 


where  the  modified  kernel  is 


(m,n) 


^PfU  asrRQ  ^  a  m6=-» 

m6=qS+Vs 


£  {(asaRV  £)(asa*V  £) 

n.rr-00  I-  r 


•  lmfi(aRqsPS)H(aRqSrS),(m)((m6+  \  qs)ebs)A(n)((v  4qs)6bs) 
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where 


/  mg\ /  m6 \ 

g6(u)  =  (asu  -  asaRqs+  ~Xv  "  aSaRqS+  ~f) 


(59) 


m6 


(|u  "  aRqSlpS)KmA('U  '  Vs^S* 


Evaluation  of  the  integrable  singularity  of  K_.  at  u=0  is  shown  in  Appendix  H. 

v>  i 


7)  Kernel  Function  K 


RD 


When  the  control  point  is  on  the  duct  and  the  loading  point  is  on  the 
rotor  the  nondimensional i zed  induced  velocity  normal  to  the  duct,  the  first 
integral  of  Eq.(3),  is  shown  in  Reference  5  to  be  equivalent  to 


V°  °  PR 


n  _  (\7)  i\7nt  ^l+a2p2 

- sin6  d0  dpD  (60) 

RD  a  p  a  a  KR  v  ' 


'7  =  LR  (PR,0RO)e  K 


R  R 


where 


KRD= 


4np  u 


NR 

-  S 

2  n=I 
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^R~0RO/aR 

VRD 
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(  3  1  B  \ 

l*R  ^  ‘  PR2 


R  R 


V§R  ‘W^W^Rn3 


'I 


dT 


RD 


RRd'  {tS+  rS  +  BR-  2rDPRC°S[6R<>-Vnt+§R-.-  *r(T  -  *0  *  5r)]  } 


On  subst i tut i ng 

i  i  00  imJ3  ®  iT, 

23  e  7  I  lm7(,k|PR)Km7(|k,rD)e  dk  (6l) 

RD  my°_®  -»  /  / 

where  5  -  eR0  “  <PD  “  «t  +  6Rn  “  aR(T  -  *D  +  SD)  . 
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and  pR  <  rD 


the  kernel  becomes 


1  y'  limit  5 _ PR  /  S _ 1  B 

RD~  ”  hrfi «  1 12  n=l  5,feR0/aR  ^  J7~J^  ^aR  '  p7  ^ 


4n2pfus 


V*D 
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I  .  7  x  J 


The  n-summation  yields 


SR  e-i<X7-m7)5Rn 

n=l 


I  'm7(|k,)PR)Km7(,klrD)e  dkdT 

-CO  /  / 


fNR  for  X^-my  =  -£jNr,  ^7— 0,^1  ,^2, 
''O  for  all  other  values 


From  the  time  t  relationship  of  Eq.(3),  =  ^7NR  =  ^8  =  ^9  =  qD 

where  qD  is  the  order  of  the  frequency  in  the  second  integral  of  Eq.(3). 
The  integral  can  be  written  as  a  single  infinite  series 


"  „  U,) 


—  ■■  (m7) 

'7  *  ,  S„  j  I  LR  7  <>VeR0>e  KRD  siW*  (64> 

y°  o  Pr 

evaluated  at  m7  -  \^-£yND  ,  and  after  the  T-i ntegration  and  taking  derivatives 
and  limits,  the  modified  kernel  is 


>7>  NRe 


-im7fD 


-iX_N„a„xr 
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i  x  0 

■  7  "\(v,VriOiS^(VV'ii,|v 
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where 


K^(2)  ‘  ■ '  ?  ts-'(2)  +  V 


and  all  linear  dimensions  within  the  braces  and  a^  are  now  fractions  of 
rotor  radius  rR^  . 

If  the  Birnbaum  modes  are  assumed  for  the  chordwise  loading,  i.e., 

(3l_)  ,  t  0-7, 1)  e  -  (X  ,n)  .  , 

V  '  i  {L»  <fR)cot  T  LR  ('>R>*i"<"-,>0«  } 

(66) 

(X?,n) 

where  L  (p  )  are  the  spanwise  loading  coefficients,  then  Eq.(o4)  can  be 

R  R 

written  as 


“  iljN.nt  (>-,,0)  (m7,n) 

'7  Vs „•  7  I  -S,  V  W  7 


*7=0 


PR  n-1 


RD  dPR 


_(m7,n) 

where  is  the  modified  kernel  after  the  9  -integration. 


(67) 
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R  RO  R 
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k  +  srVr 


(68) 


Letting  u  =  k  +  a^N^,  it  can  be  shown  that 
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•'Vd 

_(m7»n)  nr  e  ,X7°R  "'£7NRaRXD 

R°  =  '  ^Vro  6 

•  {in97(0)  +  J  [g?(u)  -  (68a) 


where 


97(“>-  ['V  +  ‘rVr  -  5]  ‘  ,U  ‘  "kVr1 


'»7<lu  -  aRV*lpX7(lu  '  ‘rVr'V 


iu(x  -a  /a  ) 


'°  'R'~R  >>((x7  -  k>J 


Since  pR  <  Rp  ,  there  is  no  singularity  in  the  original  kernel  Eq.(60) 
The  expansion  of  the  inverse  Descartes  distance  introduces  an  integrate 
Cauchy-type  singularity. 


Considering  the  L-H  side  of  Eq.(3)  in  steady  or  unsteady  case, 
certain  relations  will  exist  between  X^,  Xg,  and  Xg  and  m^,  mg,  and  m^. 
These  will  be  discussed  later. 


8)  Kernel  Function  K 


DD 


When  both  control  and  loading  points  are  on  the  duct,  the  nondimen- 
sionalized  velocity  normal  to  the  duct  at  the  control  point  is"* 


(Xg) 


i  xgnt 


*8“  ,^Ln  II  LD  ^D,PD,eD^e  KDD^XD,rD’tf,D’^D,PD’eD:^8^dSD 
SD 

*  f  .  (V  'V* 

-  S  J  f  l  °  e  Knnd6nd§ft 

Xg=°  i  2Cd  °  00  0  0 

(Xg)  Xg 

where  L0  =  duct  loading  in  lb/ft  (i.e.,  Lp° ( Cq) ) 


(69) 
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x 

RDD  =  +  rD  +  PD  "2rDPDC0S^eD  '  *V) 


The  loading  can  be  expressed  as  before  in  a  Fourier  series 

(XR)  »  (Xg.M.)  -iji0 
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Then  the  ©^-integration  involves 


2n  i(mg-p.)0D  r  2tt  for  mg  <*  (J. 


D  ^  0  for  all  other  values 


Since  Xg  =  ^ NR  *  0,  Eq.(69)  becomes 
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After  the  T-integration  and  the  successive  derivatives  with  respect 

to  Pp  and  r^,  and  nondimensional izing  the  linear  dimensions  with  respect 

to  r  ,  the  kernel  becomes 
KU 


■  im, 


'8*0 


(mo) 

^  4npfu2rR0 


{-  5  «w; 


-iSR/7N»<X0-V 
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'  2  J 
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(74) 


Examination  of  the  original  integral  reveals  that  it  is  singular 
since  R^  can  go  to  zero  when  x^  =  §D  and  Pq  =  rD  =  rq-  The  singularity  is 
of  the  Hadamard-type  (see  Reference  2)  whose  principal  value  can  be  obtained 
Furthermore,  the  expansion  scheme  for  the  reciprocal  of  R  has  introduced 
a  Cauchy-type  singularity  in  the  k-integration. 


The  peripheral  integration  over  ar*d  the  duct  chordwise  integra¬ 
tions  over  6^  and  cp^  (using  the  mode  shape  expansion  of  the  loading  LQ 
and  applying  the  generalized  lift  operator)  will  be  done  later  after  the 
last  integral  of  Eq.(3)  is  derived. 


9)  Kernel  Function 

When  the  control  point  is  on  the  duct  and  the  loading  point  is  on  the 
stator,  the  nondimens ional  induced  velocity  normal  to  the  duct  is  (cf.  Ref. 5 

for  krd) 
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2  2 

+Vs 


Vs 


sineffdeadps 


(75) 
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where 


N, 


,  ,  ’S  limit  a  a  XD~h  ,X9*'aR(T"XD+5S^'®Sn-1 
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(76) 


where  Pcn  =  -0  -  «p_  -  6C  ,  it  is  seen  that 

bU  bu  D  5n 


N|)  e~'^9+m9^Sn  „  |  NS  f°r  =  Vs’  £S  "  0,± *’  *’* 

n=l  0  otherwise 


Also  from  the  time  relationship  of  Eq.(3), 

*9  =  qD  •  (78) 

After  the  T-i ntegrat ion,  the  derivatives  and  limits  are  taken.  Then 
if  the  Birnbaum  modes  are  assumed  for  the  chordwise  loading  on  the  stator. 
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can  be  written  for  each 


as 
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where 


is  the  modified  kernel  after  the  0^-integration 


(79) 
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-  n 
K  = 
SD 
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The  expansion  of  the  inverse  Descartes  distance  introduces  an  integrable 
Cauchy-type  singularity.  There  is  no  other  singularity  since  . 

Let  u  -  k  +  aRqD 
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SOLUTION  OF  THE  SIMULTANEOUS  INTEGRAL  EQUATIONS 


1 )  Auxiliary  Analysis  of  the  Third  Equation  of  the  System 

Relating  the  three  integrals  \j,  lg,  and  |g,  for  each  value  of  X^=£ 


'7*  2  E,  e 
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The  cpp-exponential  factors  have  been  detached  from  the  kernels  and  the  remain¬ 
ders  are  designated  by  primes.  From  the  known  onset  velocities  (see  Wp  in 
an  earlier  section),  Eq.(3)  can  now  be  written  as 
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(A)  It  0 


Applying  the  cp^  integral  operator  to  both  sides  of  Equation  (82a) 
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Then  for  a  non-trivial  solution  for  all  v 
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equation 
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Comparing  this  third  ^urface_integral  equation  with  the  first  integral 
ion,  it  is  seen  that  L„  NR’^  ;s  nmited  to  the  values  l^R’^  and 


is  limited  to  L 
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(Xs,m2) 


where  “  qR  -  Xj  *  Therefore 
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From  this 


Vs  =  qR  *  0  * 

Thus  (83a)  applies  when  q  =0  and  (83b)  when  q.^0  ,  i.e.,  q  ■  X'N.  >  0 
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For  the  solution  of  Eq.(84)  each  term  of  the  infinite  v-series  is  taken 
separately  and  v  =  qR-£NR. 

By  analogy  with  I ^  and  1^,  the  chordwise  integration  over  is 
written  as 
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DD 


where 
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Now  the  generalized  lift  operators  are  applied  to  the  third  surface 
integral  equation. 
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Then  with  the  relations 
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the  kernel  becomes 
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where 
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For  the  evaluation  of  the  finite  part  of  the  integrable  singularity  of 
Eq.(86),  see  Appendix  J. 

On  applying  the  lift  operators  to  the  first  integral  on  the  R-H  side 
of  Eq. (84) 
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Equation  (84),  the  third  surface  integral  equation  for  q^O,  becomes 
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(v,m,n) 

DP 

(v,m,n) 

SD 


_(v,m,n) 

is  given  by  Eq .  (86)  ;  k  rq  is  given  by  Eq.(87), 

is  given  by  Eq.(88),  for  v  *  q  -  XN_ 

44  R  R 


and 
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When  q  =0,  tyO,  Equation  (83a)  is 
R 


_  .  (0,n)  (-£N  ,n) 

S  J  'Lp  K‘  * 


n=J  p, 


RD 


.  +  r  ,(XNR'-£V  ,(-xvHE 

%  +  I  LD  KDD  d5D 

2C„ 


(XN_ ,n)  _,(-XN  ,n)  _  (-41  ,m) 

■SfL  K  R  dp  =  W  R 

n-1  4  S  50  S  RtD 


which  becomes  (see  References  5  and  6) 

_(XNR,-XNRfn)_  (-iNR,m,n)  UNR,m) 


DD 


W 


Rtl> 


•J  K 


/«  (*£N  ,m,n)  P  (0,n)  (£NR,m,n)n 

(0'n)K_  R  +  conjugate  |_Lr  K  R 


RD 


'RD 


J/dPF 


(iNR,n)_  (-£NR,m,n) 


SD 


dPc 


(90) 


(B)  1=0 

It  can  be  shown  that  for  q  =0  and  £=0  ,  and  for  each  m  and  n 

R 

_(0,0,n)_(0,m,n)  _(0,m)  (0,n)_  (0,m,n) 

A  kdd  wd  "I  Ls  ho  dps 

ps 

(0,n)  (0,m,n) 

-J  L«  krd  %  <s’> 

PR 


hS 
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2)  Formal  Solution  of  the  System  of  integral  Equations 

The  .  integral  equations  are  solved  by  an  iterative  procedure. 

At  q  =  X'Ne,  X'=0,+l,+2,  ...  ,  for  given  order  m  of  lift  operator  mode 
K  o 

and  n  of  chordwise  loading  modes,  Equation  (1)  will  be 


WD  • 


(qR,n)  ^_(m,n) 

R  Wrr 


+  E 
X=0 


(XNr  «v.n)r  (v,m,n) 

A  LSr  (Eq.3»] 


+  S 

X=0 


„  (ANR,n)  r_(v,s,n) 

I  LS  ^PS^LKSR 
PS 


V  =  \  ~  mR  ’ 

qR  =  l'Ns  ,  X'  =  0,+l ,+2 


(92) 


Equation  (2)  will  be  at  q  =XN  ,  X=0,  +  f ,+2, . . . ,  for  given  q 

b  K  R 


w 


( XN  r , m)  (^re")  r_(v,m,n) 

s  (rs^=<?  ^  Lr  ^  >-Ki 

R  PR 


'RS 


(Eq.43, 


_(XM  ,v,n) 
+  A  R 


f_(v,m,n) 

LKd$  (Eq.52,  m5=v)J 


+ 


I 

P, 


(<0NR,n)  r  (m.n) 

LS  <p$>  LRSS  (EA.59) jdp$ 


(93) 


Equation  (3)  will  be  when  q  »*0,  whatever  X 

R 


_(XNR,v,n)^  (v,m,n)  .. 

A  Lkdo  (E,.86)]..J 


LR  LKRD  (Eq.87)]dpR 


I*  (  » n)  r_ ( v*m* n)  n 

-{  4  LKS0  (Eq.88)jdps 

PS 

(94«) 
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When  q  =0  and  X=0 
R 


_(0,0,n)_(0,m,n)  _(0,m)  (0,n) _(0,m,n) 


DD 


=  W. 


-J  S 


RD  dpR 


(0,n) _(0,m,n) 

l  Ls  *sd  dps 


(94b) 


When  qD=0  and  J&O 
R 


(^Np,-£N0,n)  (-£NB,m,n)  (£N  ,m)  (£N  ,n)_  (-£N  ,m,n) 

A  R  Rk  K  =WK-1LKK  n 

A  KDD  Vd  i  LS  sD 

1  PS 


dPc 


-H' 


(0,n)  (“£NR,m,n)  (0,n)_(£N  ,m,n)-. 

U  +  conJ  Ld  Kor^  jdp 


RD 


R  ''RD  J  "rR 

(94c) 


3)  Iteration  Procedure 

As  a  first  step,  it  is  assumed  that  rotor  and  duct  have  no  effect  on 

(*N„,R) 

the  stator  loading.  Note  that  L 


SO 


is  obtained  for  v = qR~^NR »  ^*0*1  » 


**  m^a  ~  v  # 


First  Iteration 


(0,n)  r-  \.m,n;  -,-t  r_vu,»v  -i 

)  L$0  (P$)  =  [_K$S  (Eq .59  for  X=0) J  •  [w$  (r$)(Eq.l2)J  for  all  p 


(tii.n) 


-1  r_(0,m) 


(N  ,n)  r  vm,n;  -,-i  r_V™B*"';  -) 

b|)  LS0  (ps}  "  Ks  (Eq.59  for  £=1)J  •  [W$  ( r£)  (Eq .  15)  Jfor  al  1  p 


(m,n) 


-1  r.(Np,S) 


Cj)  Then  assuming  that  the  duct  has  no  effect  on  the  rotor  loading,  L 


<V"> 

R0 


is  obtained  for  all  qR's: 


(qB,n)  r  (m,n)  -,-1  r_(qR,m) 

Lro"  <fR)-kR  C<5RMEq.2I)  j  *{LWr  <rR)<E,.4  wh«,  (6) 


PF  r  (0,n)  _(qB,m,n) 

when  qR  t  0)J  -  (Ap$)  £  |_LS0  (p$)KSR  (Eq.37  for  £=*0)  + 

P'  [Cont'd] 
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,  I  (NR’nk  ^-(qR-NR,m’n)(Eq.37  for  JM)]  } 

+  LS0  (ps,ksr  j  j 

The  loading  on. the  duct  is  obtained  in  the  presence  of  both  stator 


and  rotor. 

(0,q„,n)  _(qc,m,n) 


(Eq.86  for  v  =  qR,  L  =  0)J 


(0,m) 

r  r^n  (qR=0,Eq.l6)1  PF  (q  ,n)  _(q  ,m,n) 

{["  0  ,  ]-<»»*>  g1*. 

R 


£F  (0,n)  _(qR,m,n)  > 

’(ApS)|fLSO  (PS)KSD  (Eq.88  for  v=qR,  X=0)} 


6l)  f°r  qR  “  0 


( N  p , ”N_ ,n)  r_(-NR,m,n)  -T 

A."  *  =LKDD  (E,.  86  for  v--nr,  i.1)j 


rrr.(NR,m)  pFr  (0,n)_{-N  ,m,n) 

•{tL\D  (Eq.20  for  X=1)]-(APr)^{lro  Krd  (E 


(Eq.87,  v-  -Nr,  £-1) 


(°,n).  (NB,m,n) 


+  conJ  Llro  krd 


pF  (NR,n)  _(-N  ,m, n) 
-(Aps)SLso  (ps)ksd  ( 


(Eq.87  for  v=nr>  X-  -1)]  } 

,m,n)  n«| 

(Eq.88  for  v=-nr,  A=l)jj- 


for  qR  *  0 


(NR,qR-NR,n)  ,-_(q  -NR,m,n)  -,-1 

‘O  -  Kd  (E!i-86  for  v  ■  V  V  •"'>]  • 

r  S  (qR.n)  _(qR-NR,m,n) 

•{-<%>£  So  <s>kro  <E“-87  for  v '  W  *■’> 


% .  <V7> 


•<4ps>^Lso  Cs>Rsd 


(Eq.88  for  v  .  qR-NR,  £-1)} 


A8 
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Second  Iteration 

(0,n)  r  (m,n)  -,-1  r(0,m) 

a2)  LS1  (Ps)  =  LKSS  (Eq-59  for  X=0)J  *  {w$  ( rg) (Eq.12) 

„  pJi  (qR»n)  _(q  ,m,n) 

-  p  So  <Pr>krs  v  ■  V  fe0> 

(0,qR,n)  (qR,m,n)  , 

-  A„  K0S  <E’-50'  v  -  V  «»} 


S> 


(N„,n)  r_(m,n)  r_(NR,m) 

Lsi  ^  “  LKss  (Eq>59  for  ^1}J  “lws  (rs}  (EqJ5) 


-  <V?>  ,„<V»r.»,!) 


(pr)krs 

-  A 


(Eq .43 ,  v  =  qR-NR) 

.  (N  ,q  -N  ,n)_(q  -N  ,m,n) 

Aq  KDS  (Eq.50,  v  =  qR-NR,  £= l)j 


(qR.n)  T 

C2>  LRI  <PR>  =  LKRR  <V  (E<I-2I)] 


-1 


when  qR=0,  Eq.11  when  qR  #  0)J 


PFr  (0,n)  _(q  ,m,n) 


Wps*  [LS|'(PS*Sr"  (Eq.37  for  «) 

(Eq.37  for  jM)] 

(0,qR,n)_(qR,m,n) 


(NR,n)  _(qR-NR,m,n) 
+  LS1  <qS>KSR 


"  LAo  ’'DR 

_(NR,qR-NR,n)_(qR-NR,m,n) 


( Eq .31*  v  -  qR,  £“0) 


+  A. 


'DR 


( Eq -31,  v  -  qR-NR 


-  *”>]} 


s>  A, 


_(0,qR,n)  r_(qR,m,n)  -,-1 

A,  R  |_Kdd  (Eq .86  for  v  -  qR,  4-0)  J 


.  ru(0»"’)  -  (qR,n)  _(qR,m,n) 

•  { L 0  v°-  E<,-,6)]-  <%>  f,  Si  <s>krd  <E‘>-87-e-°> 


o  (qRf*0 


[Cont ' d3 
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pF  (O.n)  _(qR,fri,n)  «, 

-(Ws'p^S.  (PS)KSD  (Eq.88,  v  -  „R,  «»} 


e2)  for  qR  ■  0 


_(NR,-NR,n)  r.(-ND,m,n) 


LKdo  R  (Eq .86  for  v  =  -NR,  4=1)  J 

r(ND,m) 


-1 


Pf  r  (°*n)  _(~NR,m,n) 


{{  L\5  (Eq-20  for  l-')J-(Vp^iLRI  (pR)KRD 

(Eq.87,v»-NR,  4=1) 


r  (O.n)  (N  ,m,n) 

+  conj  [_Lr1  Krd  (Eq.87  for  v  =  +NR ,  4«- 1) 


PF  (Np,n)  _(-NR,m,n) 


-(aps)  pE  Lsl'R"tPs)KsD"R . (Eq.88  for  V  . -N„,  «-  1)}} 


for  qR  J4  0 


(N  ,q  -N  ,n)  r_(q  -N  ,m,n)  -| 

i  R  R  R  =  jj<^R  R  (Eq.86  for  v  =  qR-NR,  4  -  1)J 


-1 


r  pF  (<lR.n)  _(qR-NR.ni,n) 

LR1  !pR)KR0  (Eq-87  f°r  v  ■  w  W) 


■(Ap 


*  (NRPn)  „(qR-NR,m,n)  «> 

.)  E  L  R  (pjK  *  (Eq.88  for  v-q  -N  ,  4^1)j 
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LOADING  DISTRIBUTIONS 

(^NR»n) 

The  final  values  of  rotor  and  stator  loadings  LD  (rD)  and  L_  (r  ) 

R  R  5  5 

obtained  in  the  iteration  are  used  to  determine  the  blade  loading  distri¬ 
butions  of  the  respective  surfaces. 

The  chordwise  distributions  over  the  rotor  and  stator  blades  are  given 
respectively  at  each  radial  position  and  at  any  designated  frequency  by 

(a)  For  the  Rotor 

(qR)  ,r  (qR,1)  e  "max  (qR,n)  _ 

LR  <rR'V  nLLR  <rR>cotT+-£,  LR  (rR)sin(n-l)0j  (95a) 

n=z 


where  q  = any  rotor  shaf t-f reqjency  and  use  has  been  made  of  the  trigonom- 
K 

etric  transformation  x  =-0,  cos?  ,  0,  *  subtended  angle  of  the  projected 

R  DR  CX  D,n 

semichord  in  radians;  and 


(b)  For  the  Stator 


Unr) 


LS  ^S’V  ttLLS 


ir. 


)  6  nmax  (^N  ,n)  _ 

(rs)cot  —+P  Ls  R  (rs)sin(n-l)0aJ  (95b) 

n=2 


where  iNR=stator  frequencies  for  0  and  4=1;  and  x5=_®bscos®a’  ®bS  =  Su*5" 
tended  angle  of  the  projected  semichord  of  stator  in  radians. 

The  corresponding  spanwise  loading  distributions  (after  integrating 
over  the  chord)  are  given  by 

(qR)  n  <qR) 

LR  'rR^  =  I  lr  'rR’eoR^sin0adea 


*  L, 


(qBJ) 


<rR>  +  I  LR 


(qB.2) 


<rR> 


(96a) 


and 


(XN  )  n  (£n  ) 

LS  (rS)"^LS  <rS’60S)sinVea 


(*NR,1) 


(rs)  +  “  L 


(^N0»2) 
I  LS 


(rs) 


(96b) 
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HYDRODYNAMIC  FORCES  AND  MOMENTS 


A)  Rotor-Generated  Forces  and  Moments 

The  principal  components  of  the  rotor- induced  forces  and  moments  are 
listed  below.  (See  Figure  4.) 

Forces:  Fx  “  thrust  (x-di rection) 

F  and  F  =  horizontal  and  vertical  components,  respectively,  of 
^  z  the  bearing  forces 

Moments:  =  torque  about  the  x-axis 

0  and  0.2  =  bending  moments  about  the  y-  and  z-axis,  respectively 


The  elementary  forces  and  moments  of  the  various  components  can  be  de¬ 
termined  by  resolving  the  loading  force  L R  R  (rR)  acting  on  an  elementary 
radial  strip,  normal  to  the  strip,  and  taking  the  corresponding  moments 

about  any  axis.  The  forces  acting  on  a  strip  at  radius  r  of  the  ND-bladed 

I  R  R 

rotor  will  be  given  by 


AFxv  •  „5,  4  VR>  « 


cos  Dp (r)  ArR 


(R)  ^R  («1d)  iqR(Ut+®n^  R 

AFy  =  n5,  LR  (rR}  6  sJnep(r)  COS  (fit+tpRO+0n)  ArR 


/R\  ^R  (O  iqR(fit+§  )  R 
AF;R^  -  £  L-  R(r  )  e  s  i  n0R(  r)  sin  (Ot-Kp 

z  n=l  R  K  K 


RO 


where  9p(r)  is  the  geometric  pitch  angle  of  the  rotor  in 


S  i  nee 


+0  )  Ar 
n'  R 

radians. 


iqB6 
R  n 


Nr  when  qR  =  ANr,  A  *  0,1,2,... 

^  0  when  q  AN- 

R  R 


and 
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NR  (q  +1)6 
L  e  R 


NR  when  qR  =  iNR+l 

0  when  q  ^  1 

R 


Reference  1  shows  that  the  total  forces  at  frequency  iNR  acting  on  the 

ND-bladed  rotor  will  be  given  by 
R 


=  Re{l 
x  L 


iiN.nt  1  (jn.) 


NRrR0  6 


I  LR  ^rR>C0sep(rR)drR} 


,NDr..  iiN.nt  '  r  (iN  -l,n)  ,"v 

FW.R6{J^Oe  «  J£L1  «  (rR)A(n)KR) 

”  o  n«l 

(iNn+l ,n)  tZ\ .  r  . 


*  ’  (rR)A<:;)(eM)>i''eJ(rR>drR} 


N  r  n  iiN.nt  1  r  (iN  -l,n) 

R  <rR>A<  kr> 

o  n=l 


(iNR+l,n) 


<^A<")<ebR>>1"eP<‘-R>d|-R} 


The  moments  are  determined  by: 


q<R)  -  -Re  {V 


iiN.nt  l  (iN  )  R  -» 

RO  e  IL  (rR>sinep(rR)rRdrR/ 

o 


(101) 


N-rpn  iiN  nt  1  ,  r  (iN  -l,n) 

-T52  e  J{-S.  Llr  ('■R)a< ’<-ebR> 


o  rr=J 


(iN0+l,n)  -  -  (^N.-l,n)  ,-v 

+  LR  ^rR^A  (ebR0COSep^rR^+-^1i-LR  fR^  A1 

-  '■R‘N|,t';R>Al<;;)<ejR>]"<R)si"e?<Vta"6?(rR>}  rRdrR  }} 


(102) 
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and 


(d\  frNDron  'j£NDt?t  1  r  r  (JGND-l,n)  (£N.+l,n)  „ 

<)"Re{{-tr2 «  I  {-s,  L  lr  ('■R)A(n)<-ebR>-'-R  <VA<nKR> 


o  n=l 


-to  ♦£ 

(iebR,S'n6P(rR>tanep(rR>}  rR^rR  }1 


(103) 


where  A^n^(**«)  and  ( *  *  * )  are  given  in  Appendix  A. 


Thus  the  rotor-generated  transverse  forces  and  bending  moments  are 
evaluated  from  rotor  loadings  associated  with  wake  harmonics  at  fre¬ 
quencies  adjacent  to  blade  frequency,  i.e.,  at  q  =£N  ±1 ,  whereas  the 
thrust  and  torque  are  determined  by  the  loading  at  blade  frequency.  The 
steady-state  thrust  and  torque  are  determined  at  zero  frequency.  The 
corresponding  mean  transverse  forces  and  bending  moment  would  be  determined 

at  first  shaft  frequency;  in  this  case,  ^(r)  =  0  and  only  the  second 

m  R 

terms  LR  '  (r)  of  Eqs.(99),  (100),  (102),  and  (103)  are  present. 

(qR’") 

However,  in  the  case  of  the  pump-jet  system  I  is  determined 

R 

only  when  q  ==X ' N_  ,  £'  =  0,  1,  2,  ...  .  Hence,  thrust  and  torque  will 
K  b 

exist  only  at  £N  =  £'NC  (steady-state  when  l  -  V  =  0  and  vibratory 
R  b 

when  L  =  Nc,  V  =  N  )  and  transverse  forces  and  bending  moments  only  in  the 
b  K 

event  that  £'NC  =  iN  ±  1  • 
b  R 

For  example,  if  ND  “5  and  N_  =  7»  thrust  and  torque  will  exist  at 
K  b 

qR  *  *'NS  "  ^R  ec,ua^  to  ®  On  t^ie  steady  state)  and  equal  to  mNRNg  * 
integer  multiples  of  blade-crossing  frequency.  Side  forces  and  moments 
will  occur  at  t*  -  2,  i  ■  3,  so  that  ■  £NR  -  1  or  14  ■  15-1  and 
at  jI1  ■  3,  £  •  4  ,  so  that  jG'N^  ■  £Nr  +1  or  21  -  20+1. 
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B)  Stator-Generated  Forces  and  Moments 


In  a  similar  fashion  the  elementary  forces  acting  on  a  strip  at 


radius  r^  of  the  N^-bladed  stator  can  be  shown  to  be 


AF 


(c\  Js  (/n  )  i£N-(fit+e  )  s 
iS)  =  E  Lc  R  (rc)e  R  n  cosG  ( r) Ar 


n=l 


wi  th 


N, 


S  iAN.0  rN.  when  AN  *.  A'Ne 
s  e  R  n={  S  R  S 

n=l  ^ 


0  otherwise 


and 


N, 


2s 

n=l  *“ 


Nc  when  AN  = 

j  R  b 


0  otherwise 


The  total  forces  and  moments  generated  on  an  N^-bladed  stator  at  fre 


quencies  AN  will  be 
K 


i£N.nt  1  (AN.) 


Fi  )=  Re{NSrRO(rSO)e  I  LS  (rS)cOS0P(rS)drs} 


(11 


(O  fNc  i  ■  XNRf2t  J  r  (ANR-l,n)  /-»  (AN_+l,n)  / -• 

r^s)-  -*-{-4  •-Ro('-so)e  J-s,[j-s  *  <rs>A< ’Ks^s 

'  o  n=  1 


>ine|(rs) dr£  }  (II 


/,x  rNc  .  iXNRnt  1  T  (*N  -l,n)  r  (^Np+1,n)  ,-s 

tf’"  *  I  -S.ts  *  (rs)A  <-6SHS  <'s>A<> 
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The  moments  are  determined  by: 

iiN„nt  | 

o 


a‘S)  -  R'{NsrR0(rS0)ae'’eNRnt  I  Lf(rs)sineS(rs)rsdrJ 


(107) 


i£NDnt  1  r  r  (XN0-I,n) 


r rN.  iiN.nt  1  ,  r 

®  S)  -  **{{/  -V'io*  R  ft?,  [' 


o  n=l 


(rs)A(n)(-6^$) 


(^Nd+1  ,n)  -i  p  c  (^N_“l,n)  »-■> 

5  R  (^A  (6bs)]cosep|rs)^1LLs  <rs)Ai  ’ t-ebs> 


+L 


$AI"’  (9brs>  ]  ( '  ebrs)sin6p(  rs)  tan9p<  rs>}rsdrs}}  < ,o8) 


,c*  rrNc  i^N  nt  l  r_  -  (*N„+1,n) 

5<S,=  -Itofljf.  r=0(r^)P  " 


Hr-  *n.TI,n  r 

S.L4  R 

o  n=l 


(£ND+1,n)  /“\  r  c  r  (^NR"1,n)  /-\ 

-Ls  R  (rs)A(n,(eds)]co5esp(rs)+.SLLs  <r,.)A<  >(-ebs) 


+L‘1NR+(^)Ajn)(e^s)]  (iebrs)Slne9(rs)tar.eR(rs)}rsdrs}}  (109) 


where 


r1  =  — —  (nondimens ional  with  respect  to  the  rotor  radius) 


SO  r 


RO 

See  comments  on  relation  between  L  and  V  in  Section  A. 

C)  Duct  Forces  and  Moments 

The  axial  component  of  the  force  acting  on  the  duct  at  the  frequency 
£Nr  is  given  by,^ 

(*Ne)  „  (iN„) 


Dx 


R  -  r0  I  Sx  R  (*D>d*D 

2Cn 


or 


(MJ  "  (*N.) 

FDx  ro  CD  I  LDx  (xo^sineodeo 
o 


(110) 
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where  x  and  Cn  are  nondimensional ized  with  respect  to  rotor  radius  r  . 

D  D  KU 

(%) 

With  Ln  given  by.  Eq.  (97b),  the  axial  force  or  thrust  from  the  duct  is 
ux 


UV 

Dx 


=  -2r 


^  i"_ ( , 0 , 1 )  ®  _(^NR,0,n)  _  s  -i 

( 1  +cos 0ff)  +_Z)  A  sin(n-l)easineaj-decJsinQ' 

n=2 


»Lll* 


=  -2TTr  i  A 


_ ( ANr , 0 , ! )  ,  _(^NR,0,2) 


+  2  A 


} 


s  i  nor 


(HI) 


The  lateral  components  of  the  hydrodynamic  force  acting  on  the  duct 
are  derived  as  follows.  The  horizontal  (y)  component  is 


(%)  '  -27l  (%)  1 

FDy  =  roCD  COsa  J  L  I  LD  (xD’eD)sineDdeDJ5inVea 
o  o 


and  the  vertical  (z)  component  is 


(£N  )  n  r2n  (£N  ) 

FDz  =  roCD  COSa  I  L  I  LD  (xD>0D)cOSeDd0DJS,'n0A 
o  o 

(XNr) 

where  (x  ,0  )  is  given  by  Eq. (97a)  .  Since 

2TT  -iv0  r  o  for  v  ji±  1 

J  e  sin©  d0  =  -I  -in  for  v  =  +1 

o  +in  for  v  =  -1 


and 


2n 


o 


/  0 

cosO^d©.  =  ■] 

D  D  In 


for  v  jt  ±1 
for  v  =  ±1 


the  horizontal  and  vertical  components  of  the  force  become,  respectively, 


<%>. 


and 


Dy 


(£Nr) 

FDz  =  TTrocosa 


4* 


_(%,!, 1)  ]  _(ANr,1,2).  -_(Anr,-1,1)  1  _(XNR,-U). 


F„..  "  =  -Inr_coscr|  jA  "  +  *  A  "  +^A  R  Jj 


(112) 


{[a(XNr,1’°+  i  J(“r','4[a-<V,I1  i  A(ilV'',2)]} 


(113) 
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The  hydrodynamic  moments  about  the  y-  and  z-axes,  respectively,  are 


(*Nd) 


"  -  2tt  W 


'  R'  r  “1 

%  •"  ro  CDCOS*  I  L  I  L0  (xD'eD)c0seDdeDJ(VCDC0sea)sineA 


o  o 


and 


(4N  )  tt  2tt  (£N  )  *  , 

M.  =  r2  C  coso  f  I  T  L.  (x. , 6.) s i n8nd8-  l( e  -Cncos0  )sine  d6 
Dz  oD  J  L  J  D  '  D  D'  DDj  DO  cr  a  a 


o  o 


where  e-.-C..  cos0  is  the  nond imens ional i zed  moment  arm.  On  integrating,  these 
D  D  o 

become 


(*Nd) 


’R/  a  r  r-JMi  ,1.1)  .  (iN.,1,2)  ( X-NR ,-1,1)  ,  (JW-,-1,2)-, 

,0y  R  =  TTrW^A  R  +  J  A  R  +A  R  +iA  R  J 


"CdL  2 


r  ,  _C^Nr,  1,1)  ,  _(jCNr,  1 ,3)  3) 


+  tA 


+  iA 


+  f  A 


]} 


(114) 


(%)  0  r  n_(iNR,I,J)  .  _(4N  ,1,2)  1)  ,  (XN-,-1,2). 

MD2  R=-inr2c°s4DLA  R  +iA  R  -A  R  -^A  R 


r,  _(XNR,1,1)  ,  _(%,!, 3)  ,  _(1NR,-1,1)  ,  .(%,-!, 3) 


"cdLt  a 


+  7  A  -  j  A 


]} 


(115) 


When  £  =  0 


(0)  /  .  f-(°»l*l)  1  .(0,1,2)=,., 

FDy  =  Rer,TTrocosQLA  +2A  Jl 


.(0) 

Dz 


s{h 


.(0,1,1)  ,  .(0,1,2). 

2  * 


FXZ'  =  Rei+TTr.  cosoj^A  +  I  A  ]  } 


(0)  r 

mdv  =Ret 


Dy 


w  -.(0,1,1)  ,  (0,1,2)-.  r,  (0,1,1)  ,  (0,1,3)-. 

TTr-co^neDLA  +2A  J-CdL2A  +  SA  J  / 


(0)  r  .  i  r  -.(0,1,1)  ,  (0,1. 2)n  .-,(0,1,1)  .(0,1,3)-,, 

M0z  =Re(-MTrocoso}{eDLA  +  7  A  J-cJjA  +  4  A  J  } 


Note  that  the  second  index  of  A  is  equal  to  qR-^N  •  See  comments  under 
Section  A  as  to  limitations.  n 


(l  16) 
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When  i  =  1 


(N  ) 

FDy  =  -in^cosa  1j  A 


{[ 


_(NR,q-NR=J,J)  ,  _(NF?q-NR=l  ,2) 


+  2  A 


] 


L* 


_(NR,q-NR=-M)  ,  _(NR.q-NR— 1,2) 


+  2  A 


]} 


R' _ /rr(NR,q"NR  lfl^  .  1  r(NR,q"NR=1,2)‘ 


(Nr)  rr 

FDz  =  nroCOSa  lLA 


+  2A 


r _ ( Nr , q-NR=- 1,1)  ,  _(NR,q-NR— 1,2) 


+  2  A 


]} 


R'  2  |e  r-(NR»q«'NR-1'1)  1  .(Nr’1'2)  -(Nr*-1.1)  1  (NR,-l,2)n 

N„„  =  TTr_cosaq^p[A  +  j  A  +A  +jA  J 

,_(Nr» 1*1)  _ (Nr, 1 ,3)  _(Nr,-1,1)  _(Nr»“1 »3) 


"cdL7 


*1 


.  A 


]} 


<N„)  -  (X  ,1,1)  .  .(N  ,1,2)  .  .(N.,-1,1)  .  (N.,-1,2), 

C  - +7fl  -7a  -7a  j 


c. 


n  (mr,m)  ,  (nr.i.3)  ,  -<NR.-1.1)  1 J 

■^dL2  A  +  ¥  A  "  2  A  ”  4  A  Jj  '■ 

017 

When  q-NR=  1 

When  q-N  =-1 
K 

<3 


q  =Nr+1 
q  =Nr-1 
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SUMMARY 


A  theory  has  been  developed  in  treating  the  "Pump-Jet"  propulsive 
unit  comprised  of  stator,  rotor,  and  enshrouding  nozzle  by  taking  into 
account  accurate  geometry,  realistic  flow  conditions  and  hydrodynamic 
interactions  between  all  the  lifting  surfaces  of  finite  thicknesses  of 
the  system.  The  system  is  immersed  in  a  non-uniform  flow  field  of  an  ideal 
incompressible  fluid. 

The  unsteady  lifting  surface  theory  has  been  utilized  throughout  the 
analysis  and  a  numerical  solution  has  been  outlined  using  an  iteration  pro¬ 
cedure  guided  by  physical  considerations. 

Expressions  have  been  developed  for  the  various  loadings  on  the 
interacting  lifting  surfaces  and  for  the  corresponding  resulting  forces 
and  moments  evaluated  at  the  proper  frequencies. 

The  analysis  has  been  brought  to  the  point  where  the  suggested 
numerical  procedure  can  be  coded.  The  treatment  of  numerical  difficulties, 
such  as  singularities,  has  also  been  studied  and  expressions  for  their 
finite  contributions  have  been  determined  (see  Appendices  C-K) ,  This 
numerical  procedure  is  at  present  being  used  in  developing  a  computing  pro¬ 
gram  which  is  adapted  to  the  CDC-6600  or  Cyber  176  high-speed  digital  com¬ 
puter.  The  various  components  of  the  evolved  analysis  are  being  coded  for 
arbitrary  values  of  time-dependent  and  space-dependent  frequencies  and  other 
parameters  as  the  theory  indicates.  Then  by  combining  these  components  at 
the  proper  frequencies  as  the  iterative  procedure  requires,  the  correspond¬ 
ing  loading  of  all  interacting  surfaces  will  be  determined.  This  part  of  the 
synthesis  remains  to  be  completed  and  tested  for  a  realistic  pump-jet  con¬ 
figuration,  details  of  which  have  not  yet  been  provided  by  the  proper  Agency. 

Until  this  program  is  completed  and  systematic  calculations  are  made, 
no  conclusions  can  be  drawn  as  to  the  relative  merits  of  this  propulsive 
configuration  as  compared  with  a  single  screw.  Nor  can  judgment  be  made  as 
to  the  relative  importance  of  the  stator-rotor-duct  components  or  on  the 
effect  of  various  parameters,  such  as  number  of  blades,  distance  between 
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stator  and  rotor  and  their  relative  locations  with  respect  to  the  duct, 
blade  area  ratio  and  pitch  angles,  on  the  steady  state  and  vibratory 
forces  and  moments.  . 

The  present  study  is  considered  to  be  a  complete  reporting  require¬ 
ment  of  the  theoretical  analysis  of  the  pump-jet  propulsive  device.  The 
numerical  coding  when  completed  will  be  considered  as  a  supplement  of  this 
report. 
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APPENDIX  A 


EVALUATION  OF  THE  0Q-  AND  tp^-l  NTEGRALS 


f.  I(S)W  -ij  *(S)eIxC08% 

o 


where  for  m=1 


l(,)(x)  «  ^  J  (l-coscp)e'XCOScpdcp  =  Jo(x)  -  iJ,(x) 
o 

for  m  =  2 


(A— 1 ) 


I  ^  (x)  =  ^  J  (l+2coscp)e'XCOSC*:’dcp  =  JQ(x)  +  i2jj(x) 


and  for  m  >  2 


l(m)  (x)  =  ^  J  cos(m-l)cp  e'XCOSCpdcp  =  im  1  J-_j(x) 


ixcoscP^  =  :  ^ 


where 


J  (x)  is  the  Bessel  function  of  the  first  kind, 
n ' 


II.  A^(y)  =  ^  J  ®(n)e  iyCOs6sinede 


where  for  n=1 


A^(y)  =  w  I  cot  I  sin0  e” 1  ycos6d0  ■  JQ(y)  -  i J ! (y) 


and  for  n  >  1 


(y)  B  i  J  sin(n-l)6  sin©  e  'Vcos® 


d0 


(-0 


n-2 


—  LJn-2(y)  +  Jn^J 


(A— 2) 
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111.  To  evaluate 

( 

1 


-  TT  IX  COStp 

I  (x)  =  ^  l  «(m)e  cos9dcp 

o 


_  ,  TT  -iy  cose 

A|n'(y)  =  —  J  0(n)sin0cos0  e  d© 


a)  For  m=l 

l!0(x)=4[JoW  -  J2(x)]  +  ij,(x) 


for  m=2 


(2),  * 
l1  (x) 


[jo(x)  -  J2(x)]  +  !j,(x) 


and  for  m  >  2 


,;s)w .  v[-  jb(x)  ♦  js.2(x)] 


b)  For  n=l 

(1) 


A[  (y)  -  \  [j0(v)  -  J2(y)]  '  iJi(y) 


and  for  n  >  1 


.  (n) 


-  ('irt-'-[JS-3(y)  ■  Jn+l(y)] 


A,  (y) 


M 


tr 
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APPENDIX  B 

EFFECT  OF  BLADE  THICKNESS  OF  ROTOR 
ON  THE 

VELOCITY  FIELD  OF  THE  STATOR 


The  thickness  distribution  of  a  blade  section  is  represented  by  a 
source-sink  distribution  assumed  to  be  smeared  over  a  projection  of  the 
section  in  the  rotor  plane.  The  velocity  potential  due  to  the  rotor  blade 
thickness  at  a  point  (X2*rg>cPg)  on  the  stator  is  given  by 


N„  e 


R  bR 


MXS’  rs,CPS;t)Rt=  " 


n=l  - 


bR 


J 

PR 


m^r>pr,9rcP 

rrs 


1+a 


V 

RPR 


aRPR 


P-dp-dO.. 

kR  kR  RO 


(B-l) 


3f(§R,pR,e  ) 

where  M(Sn,Pn,GDr.)  =  2U - -  the  source  strength  density  determined 

R  R  RO 

in  accordance  with  the  "thin  body"  approach, 


f(§R,PR»0Rn)  =  thickness  distribution  over  one  side  of  the  blade 

R  R  RD  section  at  radial  distance  p  in  the  rotor  plane 

R 


RS 


RJ  +  rS 


xs  =  Vs  +  es  = 


E  =6  /a  =  (a  - 
RO'  R  V  R 


§Rn  “  • 


3f 


a. 


3f 


since  WD  "  eT'sTne '  W 
*R  bR  a  a 


+  PR  '  2rSPRCOS[e«0  +  »S0  ‘  ^  +  ®Rn]  Y 

(°s  ‘  ebsC0S<p«>/as  +  es  •  0  s  ‘  " 

ebR  “sV/aR  -  0  s  ec  s  n 

n~1 i 2 | 0 0 0 , Nd 

K 

,  Eq.(B-l)  can  be  reduced  to 


(•«> 


S'R. 


N„  „ 
ii  R  tt 

—  S  f  f 

2tt  n-l  J  £ 


df(pR’V 


■i 


,+aRPR 


oe 


a 


RS 


dpRd0a 
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The  thickness  distribution  f{p  ,0  )  will  be  approximated  by  a 
lenticular  section,  i.e., 

t(Pd) 


f(PR’V 


s!n  ea 


tD( pr^  a 

-2 — —  p0,  nsin20 
c  PR  bR  a 


where  t  is  maximum  thickness  in  the  projected  plane 
t 

-2-  is  ratio  of  maximum  thickness  to  chord  of  the  expanded  section 
pR0^R  is  projected  semichord 

Therefore 

2  T  (Pr>‘  PRebRSinVOS0a  ^ 

a 

The  nondimens ional  velocity  normal  to  the  blades  of  the  stator  due  to 

the  velocity  potential  (4  )  is 

h  Rt 


<ws)R  ■  •  ut^(*s)„. - 7===-  ( 


a  i  a 


's  '  Rt 


U  vl+a 


— — r  >  3x 1  r2  '  S'  r 

2  r2  S  S  ^SO  t 


(B-5) 


Substituting  Eq.(B-4)  into  (B-2) ,  and  (B-2)  into  (B-5)  and,  in  addition, 

expanding  the  reciprocal  of  the  Descartes  distance  R  as 

R5 


_L  =  I  £ 


1  ?  Im3  °n  J(x S-?R)kH. 

J<IK>me  dk 

m=-«  -eo 


(B-6) 


where  0  -  0RO  +  9$0  -  »  +  6Rp 

,  _  /  'm(«k'PR)^'k,rS) 

8" 


yields 


for  PR  <  rs 
for  PR  >  r$ 
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t  U  ''l  =  - - -  (  a  I — r  -  - —  ^ - )  T',  rf  —  (o  1  D  8.  V 1  +a“i  oz  sinG  cos  9 

\  S/Rt  rfiJTZT?  V  S  W  p  J  j  c  vkR/kR  bR  R  R  a  a 


S/Rt  r^vT^TJ 


S  “YS0  n-1  0  pR 


®  im(8B«+9(.«“nt+9r,  )® 

£  RO  TS0  Rn'  P 

m=-“  -os 


oo  j  I  J 

J(IK)me  S  R  dk  dpR  deo  (b-7) 


„lth  sR;”5»".{N« ”-**■  1-°’±,'±2" 

n=l  b  0  otherwise 


:  im0RO  "imebRcosea  im°R  ,  imCpS0  -imebSCO%  imaS 

and  e  =  e  e  and  e  =  e  e 


and  taking  the  derivatives  with  respect  to  x^  and  cp^ 


/  \  N„r_  ®  -im0,  coscp  imo  -ini2t 

(w  )  =  _ BJ _  £  e  bs  “  e  s  e 

vS/Rt  , -  m=-« 

„.,Nr 


tt  t  _ _  ima„  -im0,  cosG 

*  I  I  *T  WbR  Vl+aRPR  6  6  Sin0acos0c 

°  PR 


•  I  (iask  "  71)  <IK> 

-08  rc 


'  (xc"5R)k 

m  e  dkdpRd0^ 


(B-8) 


On  substituting  the  values  for  x'  and  §  given  before 

5  K 


iN„r„  ®  -im8,  .coscp  ima  -imfJt 

R  S  v  bS  a  5  _ 

-  Li  e  e  e 


rrsVl+a3ra  m=£NR 


•  I  I  <pR)PR6bR  Vl+aRpR  I  (aSk  '  7)(,K)m 

°  pR  -  r$ 

-i  k/a  e  cos9c  ik(a  /a$+e  )  +i(m-k/aR)oR  -i  (m-k/aR)  G^cos©^ 

•  e  e  e  e 


sin9acos0a  dk  dpR  d0ft 


(B-9) 
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The  8  -integral  involves 
a 

tt  •  m)0bRCOsea 

I  e  sin6acosVeo 


Let  u=k-a  m  in  (B -9) 

K 

n  a7  ebRCOs6Qf  i2a2 

I  e  sinVos0adeo  =  a  F(u'pR) 

°  0bR 


( B  —  1 0) 


where 


,  F  Sin(u  W  '  (U  0bR/aR)cOs(u  W  1 

F(UfpR)  =  L. - - - j 


'2NRrSaR  “  _imnt  im(1+/)aS  imaReS 
-  Li  e  e  S  e 


n2  Vl+a^r*  m=2NR 


t  -f  V  KaSu  +  V«"-7-)(IK)mF(u’PR) 

pR  bR  -®  rs 

iu^.-E+e)  _;(JL+(i+^)m)e  coscp 

e  aS  aR  S  A  aS  aS  bS  “ 
e  *e  duapD 


(B-l 1) 


where 


( 1  K) m  =  lm(iu+aRrnlpR)Km(,U+aRm,rS)  f°r  PR  <  rS 


On  applying  the  general :zed  lift  operator,  the  nondi mens ional  velocity  be¬ 
comes,  for  each  lift  operator  mode  m, 


-2aVNR  .  -imQt  lm(aS( ,+  S7>  +  aRiS) 

- Lee  S 


Wr  s(rs)  -  — Z.  7  ^ 

1  ,n=~00 
S  S  m*XN, 


’I  T  (pR)  fe:  V,+aRPR  I  (aSu  +  asV“  7)(IK>m-  F(u‘Pr) 

-00  1  C 


[Cont 'd] 
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°S  °R 

lu(S-V^ 


.  ,(m)((-  ---(?+  — )n  ebs)dudpR  (B_'2) 


On  changing  the  doubly  infinite  u-integral  to  an  integral  from  0  to  +®  , 
(B-12)  can  be  written  as 

a. 


(m) 

Vs(rs>  ■ 


■2SRrSNR 


*  -imQt 
Lt  e  e 

m=-“ 

m=£N„ 


l»(°s(l+;W«s) 

o 


|  t  «V  r;  | 

PB  bR  o 


■{lm(|uMRm|pR)Km(lu«fir.lrs)  (a$U  +  a,.aRm  -  i) 


ac  aD 

iu<4  •  if  ^ 


+  l»<l°-Vlp«)Km(lu-*R,nlrs)  (asu  -  VRm  +  pr) 


<s)((l-  -04»J 


CTS  aR 

-,u(s  -  i  ^ 


}dudpF 


(B-13) 


The  integrand  of  (B - 1 3 )  is  zero  when  u  is  zero  since  F(0,p  )  =  0. 

R 

In  the  steady  state  condition  m=.6=0  ,  the  velocity  on  the  stator  due 
to  rotor  blade  thickness  can  be  shown  to  be 


~^a «a_  r  N_  t  o  _ 

Xs  <rs>-  \  .  J  f  (P,)  ^S_v1+a»p ■> 


.  ( 0 ,  m) 

hi 


V,+as1  PR 


bR 


00  , 

•  JuF(u,pR)lo(upR)Ko(urs)-R.p.{, 


-iuA -  — +  e  ) 
e  3S  **  ,<"> 


(£  X)} 


for  »„  <  rs  . 


dudp„ 

»g  uj'  •«  R 

(B-IM 


B5 
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In  the  unsteady  case,  m  =  iND,  l  =  +l,+2,  ... 

K 


IK,,)  -'*aRrSNR 

WR  S  (rS)  =  - -  e 

1  TT2  ^/^^! 7f 

S  S 


''^rLCTs(,+  ^)+  aRes] 


*  I  T  (PR}  eT  Vl+aRPR  ”  F(u,PR)  CG8(u)-G8(-U)]dudpp 


PR  "  bR 


where 


G2^  ^  2NR(lu+aR£NRlpR)K^  (I  u+aR£NRl  rg) 


'  [asu  +  *r(v«  -  7-)] 


/i-!!.  a) 
VaS  aR  S' 


’$  -  tK) 


s  s 


(B 


B6 


-15) 
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APPENDIX  C 


Evaluation  of  Singularity  of  KDD  as  u  -*  0 

RR 


The  integral  of  Eq . (21)  is  of  the  form 


(c-1) 


where 


i^A<? 

O 


9(X)  -  lm(lX+a£Nlp)  Km(lx+a£Nlr)  B--(x)e  for  p  <  r 

B-  -(X)  =  dX+a2£N+  4)(ax+a2£N+  -^)l(m)((q-  |)ej)  A(n)  ((q-  £)0£) 

r  p 

m  =  q  +  £N 

By  L' Hospital's  rule  the  integrand  at  X  =  0  becomes 

iis  9t-2  ^  •  PfP  -  “-*> 


It  is  obvious  that 


[lrR(lX+a£Nlp)Kn(lx+a£Nlr)]  =  [ij  1  -X+a£N  I  p)Km(  I  -X+a£N  I  r)] 


and 


i  x  -  Aal 

e  X  =  0  =  e  a  X  =  0 


Xe0 


Then 


The  development  is  taken  from  Reference  2. 


Cl 
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V. 
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=  2B-  -(0){f  Km(la£Nlr)|_lm_,(la£Ntp)  +  1^,0^‘p)] 

-  L  lm(laXNIp)[K[T)_l(la£N!r)  +  K^,  ( la£NI r)]  (C-6) 

(Note  that  for  p  s  r,  p  and  r  are  interchanged  in  Eqs.  C-3,  C-4  and  C-6.) 
b)  For  x  *  0+  and  a^N  <  0 

I  ( I  x+ajWil  p)K  ( I  X+a£N  I  r )  <=  I  ((la£Nl  -\) p ) K  ( ( 1  a £N I -X)  r) 


and 


I  (I  -X+ajfcNI  p)K  (1-X+aiNlr)  =  lm(  ( I  aXN  l+x)p)Km(  (laANI+X)  r) 


m . .  r'  m'  "  '  m 

The  third  term  of  C-3  then  becomes 


lm(  (X+ 1  aZN  l)p)Km  ( (X+ 1  a  AN  l)r)] 
———  VO  . 

Therefore  Eq.  (C-3)  can  be  written  as 


-  28-  -(0) 
m,n ' 


for  n  s  r 


(C-7) 


pg(X)  9_(g.-^l"|  =  2i  ^  (IK)  I  B-  -(0) 

L  ax  Jx=0  a  mVo  m,n 


X=0 
5B-  -(X) 


+  2(.K)  I 

x=o 


a  (ik) 

±  2B-  -(0)  — - 
m,n  dX 


m 


X~° 


where  (IK)  I  is  given  in  (C-4) 


m 


X=0 


B-  - (0)  =  (a2XN+  4)  4) 1  (m)  (qeb)A(n) 


m,n 


aB-  -(X) 

— ' — — 1 1  is  given  in  (C-5) 


ax 


x=o 


aOK), 

ax 


is  given  in  (C-6) 

\  _A 


x*o 


and  the  upper  sign  is  taken  when  i  >  0  and  the  lower  sign  when  l  <  0. 


C3 


(C-8) 
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When  l  =  m  =  q  =  0,  by  the  limiting  process,  it  is  easily  shown 


that 


1 im  B-  -(0)  -  lim  ir  -  0 

ir*0  m 


m,n 


1 im  B-  -(0) 
jUO  m'n 


j^O 

S(,K)o|  ..  *2 

-  ■  ■  I  im  -*  0 

X=0  f  0  1 


1  im 

Jk-O 


(IK)o! 


X=0 


B-  -(0)  -  1 im  l  log£  -  0 

m,n  Xr.0 


38”  I  7 

lim  ( 1 K)  I  - - -  -  lim  (£+£  )log£~»  0 

£-0  °  \=0  'x=0  ir*0 

Hence  when  l  =  m  =  q  =  0 


,im  SL&Lt, S.(:A).  „  o 
\-o  x 


When 

b=» 

II 

O 

but  m  = 

q  ?  0,  it  is 

easily  shown  that 

f  1 

|  21ml 

<f) 

for  p  s  r 

O  3 

(lK)l 

=0  | 

1  1 

1  21ml 

_  1ml 

for  p  2  r 

1  im 

3('K)m 

I  - 

0 

JbO 

3\ 

x=o 

Hence  for  i  =  0,  m=q^0 


l.m  g(X)-g(-\)  _  2  Jli 
X-0  X 

1  x=o 

(™)  A  (n) 


2  jilS  <IK>m 


i 


(C-9) 


im2 


2  2 
ar  p 


(qe£)A(n)  (qej)[i  -7-7  +  am  (-7  +  -7)  ] 

r  P  r  P 

[e'l  {”>  ho;)A(S)  (=6p-ePbl  ®  (,9J)A®  <*{)]  j 


c4 


(c-10) 
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When  q  ■  0  and  Eq.  (C-9)  is  used  for  the  kernel  functions,  the  value 
of  the  integrand  at  X  =  0  is  zero  for  m  «  0.  For  m  f  0  it  can  be  easily 
shown  that  the  integrand  at  X  *  0  when  q  *  0,  in  ■  jfcN  Is 
2 

j  ^2_(a2+  4^)(a2+  4r)lm(amp)Km(amr) 
r  P 


*  {ao!  (m)  (0)A(n)  (0)  -  e£  i{m)(0)A(n)(0)  +  e£  i(m)(o)A,(n)(o)}(c-i 


I 


'  4 


C5 
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APPENDIX  D 

Evaluation  of  the  Singular  k- Integral  of  K 
The  integral  term  of  Eq.  ( 3 1 )  can  written 


_  r  F(k)  dk 
1  J  k  +  alN 

— ro 


(D-l) 


where  f(k)  -  (.k*-f-)|k|  l^dkl  rR)  [kb_,  (|  k|  ,  (|  kl  «„)] 


•  lW((m  -i)Sb)  A(n)(-kC0)e 


-ik(eD-o/a) 


This  integral  exists  in  the  sense  of  a  Cauchy  principal  value.  Therefore 


,  =  P  f[k).  :_F(-_alNl  dk  +  F(  rm  _Jk - 

1  J  k  +  aXN)  dk  ^  3  J  J  k  +  a£N 


p"  Li K)  Z-ElllM  dk 

J  k  +  a/N 


( D-2) 


where  F(-ajfN)  =  (-a2£N  +  -y*) (alN)lm(aXNrR)  {^Km_ j  (a-fNR^)  +  (a£NR0)J 

rR 

,-■>  ia£N(cp.-a/a) 

•  »(m)(q9b)  A  "  (aiNCpJe  D 

and  -F-(-a£N)  is  equivalent  to  (^)  times  the  closed  term  of  Eq-  (31)- 
For  large  |  k|  ^  |  M  j,.|  M  |  >  a/.N 

I  k|  r  .  -|k|R  .  /-k  -ik(e-o/a) 


.  (.k  ♦  .<">(-  I  k  )  A<S)  (-kC)e-ik<‘-0/*> 

r  y  rR 


••'aXN  *  throughout  Appendix  D.  The  development  is  taken  from  Reference  5. 

Dl 


t 
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The  factor 


(»k+  i)  e'|k|  (R-r>  -0.sk-- 
r 


since  R  >  r-  The  product  of  the  other  factors  also  tends  to  0  as  k 
large.  Therefore 


I  »  f  F(k)  -  F('a,,l)  dk  -  F (-a£N)  +  f  ] 

U  I.  J.  -  kl  _  u  J 


dk 


-M  k  +  a  N 


k+a£N 


S  ince 


--M 


[JWj 


dk 


k+aXN 


=  -2aXN  J* 


dk 


M  k2-aW 


,  /M  -  aXN  s 
-  '°9  H  V  a£N  > 


.  . .  F(k)  -  F(-a£N)  -  a£N, 

'  “  k  +  ajCN  dk  -  F<'aiN>'°9  (FTTS1 


The  re  fo  re 


_  (m » rn ,  n )  I 

KDR  ~  l4np,Ur 


-  imo 


fw'o  •i/l+a2rR 


•  {i"a£N(-a2£N+^)  ln(a£NrR>  [v,  (.£»«„)  *  ,  (.««,„)] 


ia£N(c0-a/«) 


lW(,VA(n)(a£NCD)[l 


J 

-M 


F(k)  -  F(a£N) 
k  +  a£N 


dk 


I 


The  singularity  In  the  k-integral 


The  Integral  in  (D-4)  can  be  rewritten  as 


,  .  (•  F'(k)  -  F'(a£N)  dk 

k  o  (k+a£N)  (k-a£N) 


whe  re 


D2 


be  come  s 


(D-3) 


(D-4) 


(0-5) 


TR-2173 


F'fk)  -  kljkr  )  [V,(k"0>  *  'W'V] 


r~\  .  ~ik(e  -a/a) 

{(ak  +  -f-) (k-alN)  I (m) ( (m  -  *)0b)  A(n) (-kC0)e 


+  (ak  -  ~)(k+aJN)  +  ^)6b)  A(n)  (kCn) 


ik(e  -cr/a) 


} 


and 


F'(aiN)  =  2a2A2N2lm(a£Hr  )  [^^(a^NRp)  +  (aXNRp)] 


•  (a^N  -  ~-)  I(m)(q9b)  A  ^  (alNC^Je 


ia£N(ep-a/a) 


At  the  singularity 


lim  f  FXk)  -  Fja/m  . 
<-a/N  l(k+a£N)  (k-a£N)J 


SF'(k) 


9k 


~  2a£N 


(D- 


k=a£N 


It  is  easily  shown  that  (D-6)  equals 


I  ^  (,eb)  A('n)(aiNCD)«iatN(€0  '  °/a) 


•  {[i(5»2«N-  -af)+  ialN(E()-f-)(a2FN-  -^)]  ^(aiNr^f-FK^a/N^)] 


+  8tNrR.(a2lN-  ~)  r(aINrR.)  [-2K^(aINR0)] 

a  ,m0(,2m-  JL)  l„(atNrs)  [V. ,  (aiNV  *  ^,<^'NRD)]} 


n 
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+  'm<aiNrR)[-2K;(aiNV] 

•  {j  (a2-£N*  ((q-2tn)6b)A(n'(-amc0)e  N<‘D  ^ 

rR 

+  iaXN(a2XN-  ^)e  (  D  I  |m)  (q6b)A(n)  (a*NCD) -CQ\  (m)  (q6b) Ajn\a£NCD)]  } 

R 

(0-7) 


slv(z) 

where  l^(z)  =  — gg  — 


skv(2) 
K^2>  =  -gr¬ 


and 


I (x)  and  A|n^(x)  are  given  in  Appendix  A. 


When  &  =  0  (m=q)  and  k  -»  0 
3F 


1  f  kl  1  (  _  _ 

3k  2aiN  =  - —rpfCm.n) 

dk  'k=ajCN  (RJ 


(0-8) 


where 


f(S,S)  -  2al®<q9b)A<S><0)  -  i  ^  [«.-§]  i®(,e.)A<;>(0) 

n 


+  i 


(n) 


■  7"  Lcd'  (q©b) ai  (°) 


-  T  I {m) (q©b>A(n) (°)] 


When  m  =  q  =  0,  SL  -  0,  k  =  0,  the  integrand  is  equal  to  zero. 


D4 
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TR-2173 


(,K)m3lu=0  =  l'"3(aRqsps)S(aRqsrR) 


for  ?S  <  NR 


\  (otherwise, 

■  VVVsN'VVr)  '  PS  and  rR  are  (E'3) 

interchanged) 


(is  ■  ‘V 


,fi(°)  =  (aSaRqS  -  3)(*rS-  7i-)A<n>((-m3t  r  0ebs)|(m,(<n’3+<>s>ebfi) 


(E-4) 


rdB-  -(u)  36-  -(u)-. 

r  m,  n  '  nyn  1 

L  du  ou  J 


&Bm  n^u^i 

m, n  I 


=  2{las(aRqS+  7^)+  ^(“VrV  72)]  A(n)(("m3+  a *  qs)  W  (  ^ (m3+qS) %) 


+("aSaRqs'  ^)(*aRqS+  TpL"1  ,im)((m3+qS)0bR)A(n)l('m3+  OO 

s 

+  i  !w  l(m)(("3+«s)ebR)A|")((-”3+^  qs)9bs)]} 


"  2{(aS  T”  -~T  -  2aSaRqs)'  (m>  ( ebRV(  n>  ((-m3  *  ^  qs)6bs) 

R 

+  i  (aSaRqS+  ^)(aRqS"  ^L"  ‘  im)  ( ^m3+qS)  0bR)A<n)CC-m3+  if  qs)Gbs) 

+  _bS  l(m)((m3+qs)ebR)Ai(n)((-m3+  af  qs)6bs)]} 


ll£f m3( I u-aR^S I  Ps)Km3( I u-aR%  1  Pr)]“  ti,m3(l  "U_aRqslPS^Km3(l"U'aRqsl  rR)]|uE< 

=  At  S(aRqSrR)Llm,-^aRqSPS)  +  S+I^rVsO  "  [Cont'd] 
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When  i=0  ,  then  the  integrand  at  u=0  becomes 


(E-8) 


When  X=0  and  m=0  ,  the  integrand  is  zero. 
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APPENDIX  F 

Evaluation  of  the  Singular  Part  of  KRg  at  u=0 

The  singularity  of  Eg. (43)  is  determined  in  a  similar  fashion  as 
before  by  performing  the  limiting  process  of  the  following  expression; 


1  im  F 

u-0  L 


9z*(u)  -  g^C-u)  i  r  Sg4(u)  Sg4(-u) 


J  L  3u 


3u  Ju=0 


g4(u)  =  'm,(N  -  VslpR)S(iu  ■  Vsl^ 


iu(e  +  a /a  -  a  /a  ) 

'  B«,n(u)  *  <f°r  V 

(F-l) 


where 


/  rn.  n  r 

BS,R<U>  -  V*SU  -  aSaRqS  -  jAV  -  aRqS  +  -p 

s  hR 

™4  *  "r  -iNR  <"1.  *  \  -  qS*  qS  '  ®V  \  *  V 

94<u>L-o”  94(-u>Iu_o 

rBg. (u)  BgJ-u)-,  ,  a  a  \r/  \  n 

Csss’ - =  2i(eS  +  Ts  -  rR  Bm,n(°) 


i  r  r  (^)  - ( _u)  -i 

+(iKm„.11[^s — =*5 — ]. 


j-8ln,4(lu-aR<lsIP|i)Kn,t(lii..Rqslrs) 

m,rPU;  L  »u 


a'",t|-u-R<1Si|;R,V<|-u-aR,,SlrS) 


■Lo  (F-2> 


I 
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)mMu=o=  'm4(aRqsPR)S(aRqsrs) 

=  ,m4(aRXNRPR)S(aRiNRrS) 


for  PR<  rS 


(F-3) 


S,5<0)  *  (aSVs  +  tX^S  '  7)  AM(H  +  ‘>S>ebR)l<m)((-V  ■£  "sXs) 


(F-4) 


9B-  n("U^ 
m,  n 


-&B-  -(u) 
m,nv  ' 


Ju  I u=0 


,-3B-  -(u)  9b-  - (u) “I  9B-  -(u) 

I  m. nv  '  _  m.  n  m  2  ^tn 

L  ou  du  Ju_q  ou 

=  2{  aS  (~aRqS+  7  )a(0)(K  +  qS)ebR),(m)((-m4  +  i^T  qs)6bs) 
R 

«r(-VrV  5  ><"’(  "  >0'<S,((  "  Xs) 


R\  S  R  S  r 


+(“aSaRqS_  7)("aRqS  +  j)l-]  if  l|m)((A+fsXs)A(n)(M)8bR) 

s  PR  s  s 

♦'  r*  |(m)((-v  S7  OO^Xh^Xr)]  } 

R  S 

'  2Kas  73  *  aR  75  -2asaRqs)'  (m)((-V  iT  qs)ebs)A(n)((m4+qs)ebR) 

R  S  5 

+  ,( W  -  £)[-  £  .  P((-V  1  ^  J 

S  R 

+  _bR  ,  Ji  qs)ebs)A{n)((m4+qs)9bR)]} 

R  S 

(F-5) 


F2 
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{  ItT  L'm,  (l“  -  VslpR)lW|u  ■  Ws1.  ’ 


-  h  [lmI)(l'','V,sl1,R)S<l'0‘aR<,slrS)]  } 


u“0 


"2  {  2  Si^R'Vs^L'm, ,-1^aRqSPR^+  'm, ,+1  ^aRqSPR^ ] 


‘  ~T~  ,m4<aRqSPR)[Km,.l(aRqSrS)+  Km,1+J  ^Vs*]  } 


for  pR  <  r$  (F-6) 


For  p  <  r  ,  p  and  r  are  interchanged  in  (F-6)  above.  (See 
R  S  R  5 


si ngulari  ty  of  K^.) 


[• 


Then  Eq. (F-2) yields  (for  l  t  0  m  ^  0) 
^94(0)  9g4(-u). 


f  - 


9o  9o  -u*0 

2i(eS  +  r  -  5r)'1.(aR1NRPR>Km<aRiNRrS) 

b  R 

•  [aSaR*V  ^IaR*R-  7]  A<n)(<™+  %>C‘bR)|(m)((-m+  V^Xs) 


rt,.<aR»R|,R>l!.(*Ri"RrS> 


+i 


{[•s  \  -  aR  \  -2asaRiNR >(B,(«“*R>0' W((*  ^Xs) 

'X  VV  tX-SV  ^  1® ((- 

S  R 

*  7T  ‘PO—UO'  ®((~  !A)»bs)]} 


-(aSW  7)(aR%"  ^^K^V^R  ),(m)((_m  +  “5p)ebs) 

S  PR  5 


[Cont'd] 


F3 


0 


i 


TR-2173 


{PRKm<VVs)L'n.-l  VNRPR>  +  'm+1  <*R*W] 


-rSlm<VVR>[lW»RVs>  *  WWs>]  }  <F-7> 


When  l  ■  0,  the  integrand  at  u=0  becomes  (for  m^O) 


/W  \ r  •  /  w  w  \ 

(rf)  {-#(*s  +  ir-ir>(n)(»V(")<-V 

ShR  5  K 


SPR 


R  S 

e 

s 


*  -fe  [♦  £  ^  <"V  '  P  (-%s>  -  TT  A‘S)  <»%>  1  <S)  <-”ebS>]  } 


and  when  m=0,  j£=0,  the  integrand  at  u=0  becomes  zero. 


1 


APPENDIX  G 


Evaluation  of  Singularity  of  K  at  u=0 


Equation  (52)  has  an  integrable  singularity  at  k=-aRq^  ' 
integral  term  of  Eq.(52)  can  be  written  as 


1  J  k+a  XN„  <! 


-  F(k)  -  F(-aR*NR) 


R  R 

(see  Appendix  D)  where 


k+aR% 


dk 


F(k)  -(.sk  -  =r  )lkl  l.dklrsJ^.dklk,,)  +KWIdklR0)] 

•'(S)((-"-fcXs><:)(-«D)«"k(V‘S'V‘s) 


(0-1) 


Fi-,AK-Vit*r7)(^%)^(n%,s,LViWDlVi(,«*«V 


*  I 


S 

t")((-"  *  ^XsX^VVl.) 


♦'•r^VV's  •  W 


and  “F(-a  XNd)  *  —  times  the  closed  term  of  Eq.(52) 
R  R  TT 


For  large  I k I  *  I M I  ,  1  Ml  >  aR£NR 
|klrs  .|klR„ 


”  (*sk*  t) 


S  *^2tt  1  k I  rs  V2|klRD/n 

'k'<Vrs> 


(S)('»r  ebs)A(n,(-kcD)' 


-,k<V*S  -os 


5  1 


■w  (-rO*1"’ <-«„>* 


!k^Ves  -  W 


which  tends  to  zero  as  k  -»®  .  Therefore, 
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S  i  nee 
r-M  “ 


Li  k«R%  ’  ‘2VNr  l  "ka-aaia»a 


(n*«> 


M  F  (k)  -  F(-a  iN.)  ,  M-apJtNR  . 

I  *  I  - rr — a--— K-  dk  -  F(-a„iNn )  log  - 5_JL) 

k+aR^NR  V  R  R;  s  \  M+a.iN.  / 


B+V \ 


Therefore 

_(m,m,n) 

KDS 


rs  eirrf7s 


4npfU  rR0 


s  s 


f  /  m  \  '  aR^R^  eD~aS  ”  as/as) 

•{•'"'rVs3.*.  *  7>  '„<Ws> 


•  0  jCN  — 

•k-l  'VNR  VS-I  <aR%V]'  <n’)((-m+  -^)ebs)A(,’)  <  VN 


r,  +  i  +  f  f<k>-F<-w  d.i 

•L1  n  9  \M+aRANR/J  _<j  k+aR^NR  dkJ 


The  integral  in  Eq.(G-3)  can  be  rewritten  as 


'k-i 


M  F'(k)-  F'(a,WJ 


_ R  R 

k2-aR2£2NR2 


dk 


where 


and 


F'(k)  =  (k  -  aR4NR)F(k)  -  (k+aR£NR) F(-k) 
F’<aR%>  ="2aRiNR  F(-.RiNR) 


At  the  singularity 


1  lm  '  F<(k)  o  5F 1  (k)  l  .  ,  , 

k-aRXNR  *■  (k+aRiNR)  (k-aRiNR)J  8k  |L„  flll  T  28R% 


k=VNR 


wi  th 


F'OO  -  klmCkrs)[Km_,(kRD)  *  K^^kR,,)] 


G2 


CD> 

(G-3) 


(G-4) 


(G-5) 


[Cont 'd] 
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“  **  * 


•  {(v  7&-  A)1' ® ((- ■  => >b>(S)  <-kcD> 


",k<VeS  -  a$/as) 


+(V+^,^>®((,+|-)9bs><")(WD)eik<e»-^  '  'S'*) 

s  s 

After  some  lengthy  manipulations,  Eq.(G-5)  becomes 
f o .  J&O  mj*0 


1  iaA(v‘s '  v“s) 


2aR%  5k 


XL'VVs’t'^V'WE  aSW  T  7+  'WvV  ijX*saR*R+  7^ 

rs  s  s 

,i;i,.ii«rs)hii«%v]|Ws'(y«"« +  7) 

s 

+lm<aR"RrS,[Km-I(aR'NRV+Km*l'aR2NRliD)]'aR'NRR0)(aSaRiNR  +  7$  )} 

♦'m^R^sX-^VVo’Hi-Z  Vr\  "t^Xs) 


•A(,’)<-aR",RCD)' 


-iaR®R(V's  -  V*sH 


+lm<aRiNRrS>[-2K;(VNRRD)]<iaRiNR)(aSaRiNR+7) 


'•  '  Wves  -  Vs) 


{ii  I  {”)((-»  +  'S“X6bs)A<"'  (8R2NRCd)”C0I  *^)ebs)A|n)(aR%l:D)} 

S  S 

(G-6) 


for  l»Q  m/O 

/  \  ni 

r$ 


(R0> 


m+1 


{[2»s + 1  7  (v‘s  -  r)]' <S)  (-™ebs>A(f;>  <°) 

1  V 


-i  4  [cD. (m)  (-mebs)A|n)  (0)  -  ^  I  jm)  (-mebs)A(n>  (0)]  }  (G-7) 


G3 
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When  m=0, 


k*0,  £= 0,  the  integrand  is  equal  to  zero. 


G4  | 

1 

1 
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APPENDIX  H 


Evaluation  of  S ingul ar i ty  of 


at  u*0 


The  singularity  of  Kgg(see  Eq.59)  at  u=0  is  evaluated  by  means  of 
L'Hospital's  rule 


1  im  q6(u)  ~  56(~u)  p96<u> 

u~*  ®  U  „  3u  3u  _  u=q 


(H-l) 


Here 


«m6(|^-aRqs|Ps)Km6(!u-aRqs|  r$)  -  Bfi-(u)  (for  p$  <  r$) 

Xn„ 

q_  =  R  m,  =  q.  +  X/-Nc  ( £,  =  0,  ±1 ,  ±2,  ...)»  and 
5  £SQ  o  b  o  s,  b 


B-  -(u) 
m,n 


(  V\ 

\  SU  “  aSaRqS  +  7J 

W  -  aSaRqS  +  ■£) 

m,+  - 

6  a 


It  is  obvious  that 


^(u)|u:0+  9*(‘u)L 


r&g6(u) 

L tst 


ag6(-u)^ 

“5H — J 


&B-  -(u) 

m  n  '  ' 


I  n 
L  ou 


&B-  -(-u)-r 


u=0 


+BS,fi<°>{  I:  [ln(lu-aFi1Nl,ll!s)K»(lu-aR"lRlrs)] 

■  ST  [ll,(|-u-aR%HJs)Km(l'u-*R1NRlrs)]  } 

u=0 

(H-2) 

<"^1^  lm6(aR£NRPs^Km6(aR£NRrS^  "  lm6(aRqSPS^Km6(aRqSrS^  for  PS  <  rS 

(H-3) 


Hi 
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•  a(Sh(»6  *  ^s)  <-*> 

o 

1  =  0,+l ,+2,  ... 

m6  =  qs  +  ^6Ns  ,  m6  =  £Nr  +  ^N$  ,  Xg  =  0)±1,±2,  ... 

SBm.n(u)  =  as(_aSaRqS+  if)1  +  if  ((m6  +  if  qsXs) 

-Si - u=0  PS  S  S 

+  as(_aSaRqS+  7)1 (m)((m6  +  if  qS)6bs)A(n)((m6  +  if  qs)®bs) 

+  (-aSaRqS  +  4)(-aSaRqS+  7) 
s  ps 

+  (-aSaRC's  +  jX'VrV  7)|W(("'S  +  IT  *sKs) 
rs  s 

s  s 

(See  Appendix  B  of  Reference  2.) 


SB-  -(-u) 
m,nv  ' 


SB-  -(u) 
m.n  ' 


Su  u=0  3u  u=0 


Hence 


SB-  -(u) 
m,nv 


Sb-  -( -u)  -1 

-S?5-2] 


SB-  -(u) 
m.n'  ' 


r 


SB-  -(-u) 
m.n  ' 


^U^i+?-2ws)'ww><">w 


[Cont'd] 


H2 
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s  Hs  s 

+ebs* (m) (zebs)A(n> (2ebs}]  (H‘5) 


where 


aR 

2  =  m,  +  —  qc 
6  aS  S 


The  first  term  of  Eq.(H-2)  is  therefore  given  by  the  product  of 
( H— 3 )  and  (H-5). 

The  second  term  of  Eq.(H-2)  is  treated  as  follows: 

Since  u=0+  and  a^q^  >  0 

'm6<lu  -  VslpS)Km6(|U  '  VslrS>  ' 

'»6((Vs  ‘  u)ps>Ki»6(<aRqs  -  u)rs> 


and 


Vslps)V(l"U’  aRqslrS^ 


'ni6((aRqS  +  u)PS)K<"6((aRqS  +  u)rS}  for  PS<rS 
The  second  term  of  Eq.(H-2)  then  becomes 

"26m,n^0^  {“  K"i^aRqSrS^Llm6-l*aRqSPS^  +  'mg+r  aRqSPS^] 

~T  l«.6(aRqSPS,LS-l(aRqSrS)  +  V,(aR<,SrS)]  }  <H'6) 

Thus  the  singularity  of  at  u=0  ,  when  ^0  mg/0  ,  is  given  by 
-3g6(u)  &g6(-u). 


L] 


u=0 


lm6<8RqSpS)V(aRqSrS,{2as(7  +  7  '2aSaRqs) 1  ^  <z9bS)A<n>  <l0bs> 

ps  s 

+  7  (7  -VRqs)(-T  -aSaRqs)L-6bSl  I"’  <26K,A<n)  (!6bS> 


[Cont • d] 


H3 
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♦  eb$|(",(26b$)Ai<;><2ebs>]} 

-2(asVs  -^)'(")Ks)A<n)(2ebs)- 

r5  PS 

1  T  Km6^aRqs'VL1m6-l  ^aRqSPS^  +  1  m^+1  ^aRqSPS^ ^ 

“  T  lm6(3RqsPS^LKm6-l^aRqSrS^  +  Km6+1  ^aRqSrS^  ) 


2  “  m6  +  qS  * 


(H-7) 


When 

1= 0  m^O  .  i 

3g6(u) 

3g6(-u)‘] 

3u 

Su  Ju= 

+ 

i  2  1  ml  [  flr 

a.  p2r2  L  bS 
b  s  s 

When 

1=0  m=0  .  ;  t 

’s  rs 


Sg6(u)  Sg6(-u) 
Su  3u 


(H-9) 


u=0 


H4 
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APPENDIX  | 

—  * 

Evaluation  of  the  Singular  k-lntegrai  of  KR[) 

The  k- integral  of  Eq.  (87)  can  be  written  as  (see  Appendix  d) 


Gfk)  -  C(-a£N) 
k  +  a£N 


dk 


(1-1) 


where  subscripts  are  omitted  and 

G(k)  =  '  2^  (ak  +  ^  1  k  1  'v{l  'K  )  LKv-1  (1  k|  RD)  +  Kv+1  (l  kl  V] 

PR 

•  A(")<<v-i)6y”N-kcD)«'k(v^ 


and 


G(-alN)  =  “  “J  ‘[closed  term  of  Eq.  (87)} 


It  can  be  shown  (following  Appendix  D )  that  for  large  (  k|  i  |  h|  >  aiN, 
G(k)  is  approximately  zero.  Therefore 


,  «  J*M  G(k)  -  G (~_aiN)  .  G(_ajtN)  ,og(  H  - 

:M  k  +  a£N  M  +  oXN 


(1-2) 


and 


_  (v.m.n) 
KRD 


Ne 


tva 


4np  U 


{+ -^  (a2£N  -  -r)l^XNpR)[Kv.|(a£NRD)  +  Kv+,(a£NRD)] 


f  RO 


PR 


.  A*"*((iN+v).9b)l*m'(3tNCD)e 


-Ui»(c  -o/a),-  .  n.aJN1 


H 


G(k)  -  G(-a£N) 


f  dk 

J  k  +  a£N 


} 


(1-3) 


-M 


The  singularity  in  the  k-integral 

The  integral  in  (| -3 )  can  be  rewritten  as 

The  development  is  taken  from  Reference  5. 


1-1 
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_  fM  G’(k)  -.G'(atN)  dk 
k  ~  J  (k+a£N)  (k-ajtN) 


(I 


whe  re 


s' 00  -  -  'v(kp,)  [v,n'  *  Kv.<kRo)] 


.  {(ak  t-^Hk-aiN)  AW((»  -  7)6  )  lW  (-kc.)e 


lk(c_-o/a) 


+  (ak  -  -^-)(k4a£N)  A(n)  ((v  +  ^)  8h)  I (m)  (kCD)e"  '  k(Va/a)} 


a'  b' 


and 


2  2 

G'(aXN)  -  -  lv(aXNpR)  [((„_,  (aXNR,,)  -  K^MNR,,)] 

-ia/N(e  -a/a) 

.  (a2lN  -  -^)A(n)((v+£N)8h)  l(m)(a£NC  )e 


At  the  singularity  k=  a£N  the  integrand  is 


lim  f  G'fk)  -  G*  (a/N)  _  aC(k)  I 

kr-aiN  l  (k+a£N)  (k-a/N)  8k  1 


-f  2aXM 


k  =  a  IN 


(I 


It  can  be  readily  shown  that  (l~5)  is  equal  to 

-  -ia£N(e  -a/a) 

l  l(m)(aXNCD)  A  R  ((v  +  XN)8b)e 


•  {  [ialU  ‘  2  “V  fa£N(€0'  f)(a£N"  J7)]lv(a£%)KJ(alNRD) 

aoB  aP  r 


+  a/NpR(a£N  -  -%)  l^alNp*)  K^(a£NRQ) 
apB 


1-2 
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+  a£NRD(3lN  -  -^  )  lv(a*Npp)  K^aJtNRp)} 


+  i  rv(«JWpR)  K'(a/NRQ) 

t  i  Ia/N(e  -a/a) 

.  {  j  (a/N  +  -^j)  (-aWCD)  A(n)((v-iN)6b)e 

aPo 


+  iaiNe 


■  ?a/N(cjj-a/a) 


(aJtH  -  “ )  [coli(m)(aJtNC0)  A^  (  (\*  *N)  6^ 


*  _b  ,  (m) 


l(m)(a£NC0)  A!(n)((v+XN)8b)]  } 


(1-6) 


where  O(z)  = 


(z) 


*;<*) 


**(*> 


a2K  (z) 


1j  (x)  and  Aa  ^  (x)  are  as  defined  in  Appendix  A- 

When  k  =  l  m  0,  it  can  be  shown  that  the  integrand  is 


n(R0) 


(V  -  - 

"TnmT  9(m»n) 


0-7) 


where 


2 
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When 


9(5.5)  -  -A(S>(veb){  ['  ♦  i  -V  (*„  -  f)]|(m)(t!>  r1 


K 


2ap 


-i  —  -^A<")(ve.)i<”)(0)  . 

n  _  .2  a  I  D 


lap 


v=0,  >£=0,  k=0  the  integrand  is  equal  to  zero. 


I  -h 


TO  W 
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APPENDIX  J 


'  Eval uation 


of  the  Singular  k-integral  of 


The  integral  term  of  Eq.  (86)  is 

,  -  -  !  r"  -iM  dk 

2  -»  k+a/N 


(J-D 


i 


where  f(k)  =  k2  [lv_,  (|  k|  RQ)  +  Iv+I  (I  M  Rj,)]  [Kv-I (1  k|  V  +  KvH  (l  kl  RD}] 

.  l(m)(-kC0)A(n)(-kC0) 

This  integral  exists  only  in  the  sense  of  a  Cauchy  principal  value* 
If  it  is  rewr i tten  as 

-a»  am 


it  can  be  shown  that 


I 


1  p"  f(k)  -  ff-aiN) 

2  J  k  +  a£N 


For  large  |  k  |  a  |  M  |,  |  M  J  >  afN 

2  e  I  k  I  R  e”  I  k|R  r"\ 

f  (k)  ^  frk2  ~  =  7~  -  I  m  (-kC)  A  n  (-kC) 

J2tt  |  k| R  J2|  k|  R/rr 

|(")(.kc)  Ats)  t.kc) 

R 


(J-2) 


(J-3) 


For  various  m  and  n  and  large  k,  the  approximate  values  of  the  I  and  A 
functions  are  tabulated  below. 

“The  development  is  taken  from  Reference  5. 

JI 


1 .  '  . 


«  - 


m  ,n 


cos(  |  k|C-  -  is!n(  |k|c-  cos(  |  k|c-  +  is  in  (  |  k  |  C  -  jj) 


cos (  i  k  Jc-  7")  +  2isin(j  k|c-  i[cos (| k|C-  ^)-  cos (1  k|  C-  J  «  0 


V 


TR-2173 


-  i  Cr"  -  f  ]  dk 


1  fH  fM  -  f(-a£N) 

2  J  k  +  a£N 
-M 


*-*[£♦  f  ]  i£l  * 


-  i  '°9<s~rsR> 


(J-4) 


where  in  the  second  term  f(k)  is  given  by  Table  ,g-2  and  where 
f(-a£N)  =  a2l2H2  [l^,  (aXNRD)  +  lv+,  (**!%)] 

.  [kv_ j  (aiffRp)  +  K^,  (a£NR0)]  l(m)(a£NC0)  A(n)(a*NC0) 


In  the  case  n  >  i  where  f(k)  -*  0  as  k  —  ®,  the  kernel  becomes 


_  (v.m.n) 


DO 


<mpfUro 


{j  *W['v-l  (8iNV+'v.l  (a“R0)][Kv-l  (’""V^v+l  (a“% 


l(")(a£HC„)A('n,(.iNC[))  [l  +  +  log  (^ffj)] 


M 


j.  r  dk  \ 

2  J  k  +  a  N  aK  J 


When  n  =  I,  there  are  additional  terms  which  may  involve. 


const,  {j  ♦  f  ]  -  const. 

“OB  n 


(J-5) 


and 


J3 
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[r\f  “I  9<k>dk  /  r“  g(-kldk  r' 

■‘•U  +J  J  i?aiN-=  const*  l  -  J  rrtis  +  J. 


M  k  ~  aiN  JM 


Qfkldk 

k+a£N 


where  g(-k)  =  _ . 


~ i2kC  -  .. 

+e  ,  m  odd 


-  i2kC 

-e  ,  m  even 


•2kC  -  .. 

f  -e  ,  m  odd 

9(k)  ■  {  I2kC  - 

+e  ,  m  even 


The  integrals  are  evaluated  below.  By  means  of  the  substitution 
X  =  2C(k+ajfcN) 


r“  e‘2Ckdk  =  -i2Ca£N  p  e*XdX 

M  k  +  aXN  2C  (M+aJIN)  X 


-?2Ca£N 


Ci  |" 2C(M+aiN)]  -  isi  [2C(M+aiN)]  } 


where  C?(x)  =  -  f  — 
si(x)  -  -  f  ^ 


cos  XdX 

S  i  nx 

X  ~ 

x 

s  i n  XdX 

-cosx 

X  -  ~ 

X 

for  x  »  I 

for  x  ».  J 


(See  Jahnke  and  Emde:  Tables  of  Functions.  Dover  Publications,  New  York, 

19450 


Therefore 

r”  e'2Ckdk  -i2Ca£N  (  -sin  [2C(tt+aXN)"|  +  icos  r2C(H+aXN)l  ) 

H  ( - 2C{LiU)  - - j 


.  i2CM 
1  e _ 

2C(M+aXN) 


(J-6) 
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f'(aXN)  -  -Ui(alH)3  i;(aXNRD)KJ(a£NRD)!(fn)  (a£NCD)  A(n)  (aXHCp) 
It  can  be  easily  proved  that 


dJUIkl  ^  2a£N  = 


3k 


k=a£N 


-2i{2a£N|^(a£NRD)K^(a£NRD)+a2£2N2R0  ^'(aXNR^K^aANRp) 

+  a2/2N2RDl^(aXNRD)K”(aANRD)}  {l (m) (a£NCD)A(n)  (alNCp)} 


+  i 


iaXN  {l^(a^NRD)K^(aANR0)} 

.  {l(i;,)(-a£NCD)A(n)(-aXNCD)  -  I  (m)  (a£NCD)A(n)  (aXNCp) 

+  i2a£NC0[l(i;,)(aXNCD)A1(f;)(aXNCD)  -  I  fm)  (a£NC[))  A(n)  (aANC,,)]  } 

(J-9) 


When  1  =  0 


,  x  _2vCn 

|iM  2atN  -  — ^ 

3k 

k=a£N 


[l{m)(0)A(n)(0)  -  l{m)(0)A1(n)(0)]  (J.10) 


where  (see  Appendix  A) 


1 

”  2 

for 

•  m  = 

«  1 

(m) 

(0)  -  \ 

.  1 

.  0 

for 

for 

m  = 

m  > 

1  ,2 

2 

■  (m) 
'l 

(0)  = 

1 

for 

m  = 

2 

0 

for 

m  > 

2 

1 

for 

n  = 

1 

1 

t  1 

I  5 

for 

n  = 

1 

(n} 

(0)  =■ 

1 

5 

for 

n  = 

2 

A<",) 

(0)  = 

]  1 

1  4 

for 

n  * 

3 

0  for  n  >  2 


0  for  all  other  n 


J6 
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APPENDIX  K 

Evaluation  of  the  Singular  k-integral  of  K, 
The  k-integral  of  Eq.(88)  can  be  written  as 


SD 


»  G(k)  -  G(-a  JN  ) 

i  *  J-  R  R 


dk 


k+aRXNR 

M  G(k)  -  G(-aR£NR) 

“  k+a  £Nd  ■>  -i  i  •  “p~”D' 

-M  R  R  R  R 


dk  -  G(-  ^r)1o9(^-JO) 


(K-l) 

(K-2) 


where 


G<k>  ’  -  isr  (%k  -  7-)lkl  'vdkiPs)rLKv-i<ikiRD> +  W'lV] 

S  HS 

.A(S)((.v.L.)ebs)l(")(.s).ik<V's-VBs) 

s 


and  G(-aR£.NR)  ■  -  ^  {closed  term  of  Eq.(88)j  . 

Therefore 

-i  vcts 

k  (V>5'") «  V _ r ,  !»  %  (, ,  u  „  .v-'A'Vs  -  V*s> 

50  ',n'fuSrR0  *  *s  2  ^ s  R  R  v 


'  Iv(Sr'Vs)[Kv-I(VNRRD>  +  WWd*] 


['  ♦  b  -  ] 
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M  G(k)  -  G(-aRXNR) 


k+VNR 


dk 


(K-3) 


Kl 


where 


G'<ki  -  -  2^  lv<k<’S>[Ku-l<kRD>  +  Kv+l<kRD’] 

•{(ask-7)(k-aA)A<S)((-v-^Xs)'(S)(-S> 

ps 

+(ask+7)(k+aRivA(n)((*v+yebs)|(m)(kcD)e 


!k<  VeS  ’  °S/aS) 


■ik(eD_eS  '  W 


a?X2Na 


G'<W=  -~sir-  ^VVs^WVW  +  WWtf] 

r/  \  r\t/  a„iNBN.  N  ,-v  "iaDiNt>(eB"ec  "  CTc/ac 

7-Vn)((-u  * -^ks)1  (aRJVo)e 


11  m  f  G'(k)  "  G'(VN^  l  _  3G*lii 


k-aR£NR  L  (k+aR£NR) (k-aRiNR) 


.  -  2  a 

k  UA  R  R 


(K-5) 


Equation  (K-5)  is  equal  to 

i  t~\  r~\  //  aD£N„\  \  - i  a_  ^N_  ( en_ec  -  <?  /a  ) 

VuV1  ((-  +  V^ 

'{[  I  aR*V  I  a  ca  "  '  *R%(W  aj  )(aRfllR+  a  p>)] 

s  s  s  s 
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+  i  'vfWsKVW 


+  iiW 


•'WVS  -°S/aS) 


(W^K'^'Wo’ 


-<"')((- ^K) 


R  R  a  os/L  D 
SKS 


|W(aR«RCD)A|n>((-v 


J]  } 
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where 


km 


*km 


km 
KM  ■  -fe- 

a2Kv(2) 

Kv<z)  c 


I (x)  and  a|"^(x)  are  given  in  Appendix  A. 


V/hen  kgl°0  ,  it  can  be  shown  that  the  integrand  is 


I  VV 
"  (V"' 


g'(m,n) 


where  g'(m,n) 


>)(-vebs){[i  -  i  -V  (V«s  -  J)]  '(S)«» 


2  Vi 
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or 


g'(m,n)  - 

i(m^(o)A(n)(-vebs)  [- 


1  + 


2a  -2 


'SPS 


(vs  -  $1 


e. 


-  t  J7T2LcDlS">'0)A(,’,<-v6bs)  '  if  |(m)(0>Ain) 

zasps  s 


When  , £*=0 ,  the  integrand  is  zero. 
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APPENDIX  L 


EFFECT  OF  RACE  OF  STATOR  (S)  ON  ROTOR  (R) 

In  the  course  of  a  re-examination  of  the  theoretical  development  of 
Reference  9,  it  was  found  that  the  behavior  of  the  velocity  field  for 
points  inside  the  propeller  race  is  quite  different  from  that  at  any  other 
point  in  the  field  around  the  propeller.  The  existing  theory  and  program 
dealing  with  the  propeller-induced  velocity  field  have  therefore  been  modi¬ 
fied  to  include  the  region  of  the  propeller  race.  The  wake  effect  of  the 
stator  on  the  rotor,  designated  by  AW  /U,  has  been  developed  and  incorporated 
in  the  present  program  to  be  used  whenever  there  is  no  available  wake  survey 
at  the  rotor  plane  in  the  presence  of  the  hull  and  stator. 

The  WR  induced  velocity  at  points  on  the  right-handed  after  rotor 
by  the  presence  of  a  "left-handed"  forward  stator  is  given  by 


W. 


a) 


u'WV')"  |  I  x504ps  (WV 

s  s 

a  >  ~i*[aD(T'-xD)-er_] 


+  ixnRt 
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nR  -® 


R'  Sn' 
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where  X  -  £NR,  X  =  0,1,2 . 

A 
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^RO  aR  ^*R 
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R  3  R 
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ln  RSR'{(T,'5S>a+rR+PS'2r«PSCOS[6S0'KPR0'  QRt+SSn]/  ' 

Let  t'-x_=t,  and  0  = -Q  t  .  The  T-integral  then  yields 
R  R 


O  t  Jk  \  /  |  \ 

'T'i  e  k  *h  *  pT-  5vXR^)psdpsd5s<iT 


where 


R5R=  {(T+XR-?S>S  +rR+PS-2rRPs“s[eS0^R0  +  6  +  5Sn3}' 


Then 


a,T  °P  *IX(aRT-eSn) 


3x 


■I 


R  - 


^  4^)(^)psdpsdVT 


But 

Therefore 

3't 


da 

K 


=  -I 


o  +iX(aT-§Sn)  f  a2  }  S2 


(a  l!_  +  I - 1 - )(¥i-)psdp  d|  dT  (L-2) 

'S&x|  Pg  3x^30Sq^^sr  S  S  S 


Furthermore  for  points  inside  the  propeller  race,  Laplace's  equation 
written  in  cylindrical  coordinates  takes  the  form 

S  (IX  h  I?  (P  I?  D+  ?  iF  ©-  -  r  6(™«is)6  VPS),(  VW®*V 

R  o  *> 

where  6(  )  is  the  Dirac  delta  function. 

Thus  whenever  the  field  point  coincides  with  the  helices  of  the  wake, 

S'  (s)=-T  6(T+*R-?S)6<rR-p5)S(6S0*%>  *  ®  *  5Sn>-L^I^PS  t^X*) 


s  so 


rs 
(1-3) 
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The  induction  can  be  expressed  in  a  Fourier  series  expansion  as 

Au(')  „  TT  n 

1  r  R  -in© , 


■jr-  -  E  CV 
U  nc"°°  ** 


Let  0=6'  .  Then 
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=  n5-  — |?4  e 


sin6adeo 


Since 


NS  'n®5n  r  Ns  w^en  n  *  ^  Ng,  ■  0,  ±1,  ±2 . 

E  e  ■  t 

n=1  v  0  otherwise 

Equation  (L-4)  becomes 
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and  assuming 

(^•)  ,  Q  »  1  y'  (X,n)  - 

LS  (rR*0a)  "  Wfj?)  LS  (rR)©(") 

where  0(n)  represents  the  Birnbaum  chordwise  modes,  then  the  integral  part 
of  (L-5)  becomes 

i\  —  a 

maxn  iXaDe  JnaD  -in6.Bcoscf>  i  (X+n)ac  a  S 

I  -  E  S  LiX>n)(rJe  R  Se  Re  bR  “e  Se  S 

n-1  \=0  S  R 


•  A(5,(("^0*^)>bs) 


(L-6) 


Taking  the  lift  operator  at  each  m-order  and  nondimens ional i zing  with 
respect  to  rQ(rotor  radius),  Eq.(L-5)  can  be  expressed  as 
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(L-7) 


where  n  *  £jNs,  £j  *  0,  ±1,  ±2,  ...  . 

X  -  £Nr  ,  £  ■  0,  +1,  +2,  ...  . 

(it  can  be  shown  by  a  similar  approach  that  the  second  term  on  the  right- 
hand  side  of  Eq . (L-3)  does  not  contribute  to  d3/dx^(~).) 

In  the  steady-state  condition,  £j*0,  and  retaining  only  the  tr*Q  and  1 
terms  (i.e.,  X“0  and  X«NR)),  Eq.(L-7)  becomes 
(0,m) 
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APPENDIX  M 

THE  VISCOUS  WAKE  OF  THE  STATOR 

In  a  pump-jet  propulsive  system  the  rotor,  being  located  in  the  race 
(wake)  of  the  stator,  operates  in  a  real  fluid  and  hence  should  include 
both  the  potential  and  viscous  effects.  In  the  absence  of  wake  measure¬ 
ments  in  the  plane  of  the  rotor  when  the  stator  is  in  place,  it  is  necessary 
to  take  this  into  account  theoretically.  The  potential  contribution  has 
already  been  dealt  with  in  Appendix  L.  The  effect  of  the  viscous  wake  is 
approximately  considered  by  the  Kemp-Sears  method  described  in  Reference  10. 

The  configuration  of  viscous  wakes  of  propeller  blades  is  approxi¬ 
mated  from  single  airfoil  experiments.  The  unsteady  force-and-moment  on  a 
downstream  blade  passing  through  such  wakes  is  then  calculated  on  the  basis 
of  the  theory  of  isolated  thin  airfoil  in  nonuniform  flow.  The  same  approach 
has  been  adapted  to  the  unsteady  lifting  surface  theory. 

Silverstein,  Katzoff,  and  Bullivant,"  have  shown  that  the  half-width 
of  the  wake,  Y  ,  may  be  calculated  from  the  following  formula 

Y  =  0.68  CD5  c(x/c  -  0.7)*  (M-l) 

where 

c  *  airfoil  half-chord 

x  *  distance  measured  along  the  wake  axis  (free-stream  direction) 
rearward  from  the  center  of  the  airfoil 

Cp  ■  the  airfoil  profile-drag  coefficient 

1 2 

NOTE:  Cp  will  be  calculated  according  to  Hoerner's  method. 

For  convenience,  a  new  coordinate  x“  along  the  wake  axis  is  Intro¬ 
duced  in  Eq. (M-l)  : 

x"  ■  x  -  0.7c  (M-2) 

f  0  *y 

Kemp  and  Sears  have  shown  that  in  terms  of  x  the  wake  half-width  and 

Ml 
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the  velocity  at  the  center  become 
Y  ■  0.68  tjl  c(Cdx'7c)j 
u  c/V  =  -(2.i*2C0^)/(x*/c  +  0.3) 


(m-3) 

(M-4) 


and  that  the  velocity  profile  to  be  used  is 


c 


exp 


(M-5) 


Since  the  propeller  blade  moves  along  a  line  oblique  to  the  x 
(or  x")  axis,  it  is  convenient  to  introduce  oblique  coordinates  x',y'  as 
shown  in  FigureS.  The  relation  between  x",  y  and  x'.y1  is  given  by 

x  =  x'  -  y'cosQ^  ,  y  =  y'sin8^  (m-6) 

P  P 

(The  superscripts  S  and  R  refer  to  stator  and  rotor  blades,  respec¬ 
tively.) 

Since  the  wake  is  narrow  in  the  region  of  interest,  see  Figure  5» 
y'/x1  is  small  in  the  wake  itself,  and  one  may  write,  approximately. 


RJ  x1 


^sin6S 


(M-7) 


Then  the  wake  half-width  and  centerline  velocity  are  as  follows: 


•■«K(£x£)] 


i 


( M-8) 


u 

c 


where  c^ 
The  veloci 


-  °-3) 

is  the  total  chord  length  of  the  stator, 
ty  profile  from  Eq.(M-5)  is  now 
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r(0  -  f-  -  Y) 


r  /  r  s  i  nB^v  2  /  \  2-« 

5"  *  exp[‘T{ — y  ^ )  (e  "  T  ■  Y)  ] 


( M— 11) 


where 

0  —  angular  coordinate  of  the  stator 
Y  —  angular  coordinate  of  the  rotor 
r  —  radial  posi tion 

Equation  (M-l 1 )  can  be  expanded  in  a  Fourier  series  in  terms  of  (6-Y) 


—  =Z/(a  cosn(0-Y)  +  b  sinn(0-Y)) 
u  n  \  n  '  '  n  / 


(M-12) 


or 


where 


—  =  S  (a  cosncp  +  b  sinncp) 

u£  n  n  n 


9  -  0  -  Y 

2tt 


an  -  ^  J  (H_)cosncpd9  (NR=  no.  of  blades  of  rotor) 


(M-13) 


(M-l  4) 


(M-15) 


o  c 


N  2tt 

bn  =  2^  I  (^)sinncpd<p 
o 


(M  - 1 6) 


The  velocity,  u£  ,  is  in  the  direction  of  x  ,  which  makes  an  angle 

(0^+0^)  with  the  after  propeller  blade  so  that  the  component  giving  upwash 
P  P 

at  the  blade  is 


^  •  i-  Sln(es+e») 

u  vS  u  P  P 


(M- 1 7) 


and  since 

£ 
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s  i  n0* 


then  from  Eqs.(M-9)  and  (M-17), 
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(M-18) 


where 


x*  cR  (e  p.S  ■  xR  VS \  n  , 

$  S  V.cR  P  c& 


Choose  — r  =  0,  which  means  the  point 
c 


rotor  blade.  Then 


(£) 


CSC  0°  -  0.7 
p 


(M  *  1 9) 


is  at  the  mid-chord  of  the 


(M-20) 


The  viscous  wake,  then,  can  be  expressed  in  the  following  form; 


.<«!> 


c 

U 


(a  cos  ncp  +  b  sin  n  cp) 
'  n  T  n  T/ 


(H-21) 


where 

q  =  2n  (M-22) 

cp  =  0  -  Y  «  20  (M-23) 

The  left-hand  side  due  to  unsteady  wake  in  the  PPEXACT  (Propeller- 
propeller  Exact)  program  (Reference  1 )  i s ,  in  lift  operator  form, 

-fa*"')  (q)  .  r  ,-x 

-X  (r)-^(r)  e-1"0  lW(<)  (H-2k) 

where 

1  77  iq0.rcos^ 

|(m)(qebr)  =  W  J  e  (M-25) 

o 

8(1)  *  1  -  cos9a 

§(2)  =  1  +  2coscpa 

8(m)  =  cos(m-1)cpff  for  m  >  2 


Thus,  the  resulting  unsteady  force  and  moment  or  unsteady  side  force  and 
moment,  at  the  specified  blade  frequency,  can  be  determined  as  in  the 
PPEXACT  program.  These  viscous  effects  are  then  superposed  on  the  results 
from  the  potential  flow. 
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